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ABSTRACT 

Information Theory Inference, better known as the Maximum Entropy Method, is used to 
infer the lateral density distribution Inside the earth. The approach assumes that the earth consists 
of indistinguishable Maxwell-Boltzmann particlep populating infinitesimal volume elements, and 
follows the standard methods of statistical mechanics (maximizing the entropy function). The 
GEM I OB spherical harmonic gravity field coefficients, complete to degree and order 36, are used 
as constraints on the lateral density distribution. The spherically symmetric part of the density 
distribution is assumed to be known. The lateral density variation is assumed to be small compared 
to the spherically symmetric part. The resulting information theory density distribution for the 
cases of no crust removed, 30 km of compensated crust removed, and 30 km of uncompensated 
crust removed all give broad density anomalies extending deep into the mantle, but with the density 
contrasts being the greatest towards the surface (typically ±0.004 g cm"^ in the first two cases and 
±0.04 g cm‘^ in the third). None of the density distributions resemble classical organized convec- 
tion cells. The information theory approach may have use in choosing Standard Earth Models, 
but, the inclusion of seismic data into the approach appears difficult. 
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INFORMATION TIIKORY lATERAL DKNSITY 
DISTRIBUTION FOR FARTII INFFRRm) 
FROM (iLOBALGRAVITV FIFU) 


INTRODUriTON 

The problem luUIrcsseii Fere is inferring via information theory the lateral density variation 
inside the earth from tlie obscevetl external anomalous gravity Held, Information theory is \ised 
beeause It Is the least subjective way to deal with inverse problems (Balerleliu 1971). The motiva- 
tion for this study is the relationship of the lateral density variation to tectonics and convection. 

The natiuv of the problem is the following. The lateral density variation inside the earth 
generates the observed gravity anomalies. Hence information about the lateral density variation 
is provided by examining the anomaHes. However, the observed gravity anomalies cannot be 
inverted to recover tl\e actual density variation. The pioblem is nommique: there am an infinite 
number of density distrlbutioivs which can generate the observed gravity field. This is unfortunate, 
since it is desimblc to know the lateral density variation, especially with regard as to how it relates 
to tectonics and convection (Phillips and Lambeck, 1980). 

The nonuniquenoss can be dealt with by various approaches in orilcr to obtain insight into the 
physics of the earth. Of these modeling Is by far the most common approach. Hem extra assump- 
tions are introduced until the solutlott to the problem becomes unique. Kaula (19()3), for instance, 
assvimes that the shear sintin energy of the mantle is minimised. This key assumption, and other 
minor ones, together with the constraints of the observed gravity field determine a unique latend 
density variation, Phillips and Uunbeek (1981) review many papers which use modeling. Another 
approach is the Backus-Gilbert meiliod (Backus and Gilbert, 1967. 19{)8; Parker, 1977), which 
studies all possible solutions consistent with the given data, This study is called the geophysical 
inverse problem (Backus and Gilbert. 1967, p, 249), The Backus-Gilbert method has been used 
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extensively in seismology (e.g., Jordan and Franklin, 1971). Burkliard and Jackson (1976) apply 
the metiiod to gravity data. There are also other approuches to data inversiom see, c.g„ Parker 
(1977), Sabatier ( 1977), and references contained therein, 

The approach to the lateral density variation used here is that of information theory (Jaynes, 
1957, 1963, 1967). This approach is commonly known as the Maximum Entropy Methotl, or 
MEM for short, A better name would be Information Theory Inference (ITI for short), since 
information theory is its basis and the name avoids conihsion with thermodynamic entropy 
(Baierlein, 1971, pp. 47 3-478). It will be called ITI here. 

ITI is a probabilistic approach to nonuniquencss. Each possible answer (labelled i) to a non- 
vmique ])roblem is assigned a probability Pj that it Is the coriect answer. The probabilities Pj are 
assigned numerical values so as to maximize Shannon’s (1948) Information measure 

MI«~SPjlnP, (1) 

subject to the constraints of the known data, Ml in (I) stands for ‘'Missing Information,” l.c., 
the amount of information needed to determine which answer is correct (Baierlein, 1971, p. 64). 
In practice the expectation vahie of the desired unknown quantity is taken as tlie infericd answer 
to the problem. T he information tlieory approaclv thus provides a solution to what wiU be called 
the geophysical inference problem: picking one answer out of a number of possible answci's as the 
most likely to be true. ITI may hence be regarded as a complementary method to the Backus- 
CUbert method (Gull and Danlell, 1978), which investigates the geophysical inverse problem. 

The rationale for using ITI is that it picks the “best” answer, dictated by the information at hand, 
out of the many possible answers. "Best” here means “least subjective,” i.e. the number of 
unconscious assumptions in choosing an answer is minimized (Tribus and Rossi, 1973), Sec 
Baierlein (1971, pp. 11-89) for an excellent introduction to ITI. 


IT! (or MEM) )«as bct'u used svith great success In several different fields. One is statistical 


mechanics (Jaynes^ l9S7;Tribus, 1961; Kat?., I967;and Balerlein, 1971). In fact, (1) is the 
iMitropy function of statistical mccimnics. The only difference between ITl and statistical mech- 
anics lies in ITEs powerful information theory foundation, which allov/s the approach to be 
applied to a wide variety of problems, and not just to statistical mechanics. It has been so applied 
to spectral analysis (Burg, 1967, 1968, 1972) and to radio brightness maps of the sky (Gull and 
Daniell, 1978). In solid earth geophysics ITl has been applied to the spectral analysis of polar 
motion (e.g„ Smylie et al., 1973; Graber, 1976), and to the mdlal density distribution of the earth 
(Rietsch. 1977; Rubincam, 1978, 1979; Graber, 1977) and of the planets (Koyama, 1979). 


The present work is an extension of Rubincam (1979) to the lateral density variation, ITl is 
used to Infer the lateral density structure based on the spherical harmonic coefficients of the 
gravity field and on an assumed spherically symmetric density distribution (also called the radial 
density distribution). The gravity field coefficients are tliose of GEM lOB (Lerch et al., 1981), 
complete to degree and order 36. The hydrostatic equilibrium bulge is subtracted out of the 
ft » 2,4, m * 0 terms using the hydrostatic coefficients of Nakiboglu (1979). Also subtracted from 
the GEM lOB terms whei; the need arises are tlie gravity field coefficients of a crustal model. 

Two different crustal models are used: one a 30 km-thick isostatically compensated criist and the 
other an isostatically uncom;)ensated crust, also 30 km thick, Carl Wagner supplied tlie spherical 
harmonic coefficients for tliesc models (Wagner, private communication, 1976). Both sets of 
coefficients arc complete to degree and onler 36. The radial density distribution is the average 
structure Parametric Earth Model (PEM) of Dziewonski et al.,(1975). Further, the earth is assumed 
to be a sphere and that the Nteral density variation is small compared to the radial density distribu- 
tion, 

The principal results are as follows. The intbrmation theory density distribution can be 
written as a spherical harmonic expansion. The equation for the density variation is similar in form 
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to that or the equation whicli gives density contrasts due to lateral temperature diffemnees. 

For the cases where no crust or the 30 km thick compensated crust Is removed the density con- 
trasts are greatest near the earth*s surface and have typical magnitudes ofiO.004 g cnf^ . The 
density contrasts are also greatest near the surface for the case of 30 km of uncompensated crust 
removed and are typically a factor of 10 larger than in the two other cases, In all three cases the 
density contrasts decrease witir depth but slgnlUeajit density anomalies still extend deep into the 
mantle, None of the three density distributions look like classical convection patterns, i.e„ 
oi^ganized cells with columns, of low density where material Is rising and columns of high density 
where material is sinking. No attempt has been made in any of the cases to compute stresses or 
stress-differences. 


DERIVATION OF THE INFORMATION THEORY 
DENSITY DISTRIBUTION 

The information theory density distribution is derived from the following considerations. 

It is l^t assumed that the data consist of the known values F^j for Q integrals of the form 

Fq s J pir) f^ (?) dv, (q « 1 , 2 , . . . , Q), (2) 


where p(r) is the earth’s density distribution, fjj(r) is a function which depends on position ? 
inside the earth, dv is a volume element, and V is the volume of the earth, Examples of integrals 
of this form are the mass and moment of inertia of the earth. The gravity field coefficients also 
have tius form, since they can be written (Plullips and Lambcck, 1580, p. 30) 


i 


yp (r) r^ Yfi,„i(Q.\)dv 
(2fi+l)Mjj a| 


fiini 


(3) 


where and are the normalized coefficients of degree fi and order m, 

r = I ?1 , and the Yg„,j(0,X) are surface spherical harmonics using Kaula’s (1967) 4»r normalization, 
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wUh 0 being colutilude ami X longiluclc, Mj; aiul a|« are the mass ami radius of the earth, respec- 
tively. Obviously in this case * Cg,„j and 


n?.- 


(4) 


The next eonsklcration in using ITI is to set up the earth models which constittUc the various 
possible answers to the problem. The task Is then to choose the “best*’ model based on the data 
of the form (2). The earth models are set up as follows: the earth Is divided up into infinitesimal 
cubes, all with equal volume dv « dxdydz. Bach cube is labelled with the running subscript j. 

The position vccto-r from the center of the enrih to the jtli cube is^, The cubes are populated with 
Indistinguishable particles of mass m. Bach earth motlel i has Ujj particles in the jth cube where 
Ujj is an integer >0. The density of the earth at position in the i th model is then « njjin/dv, 
The integrals (2) for each earth model i become the sums 


f 


qi 


/njiiuX ^ ^ 

f \lk”/ f 


(5) 


The expectation values 


P > a V n p , 

j I 


(6) 


are assumed to constitute the observed values of F^, where P| is the probability that tlie ith model 
is in fact the correct model. 


The problem so formulated is analogous to the standard statistical mcclianics problem of 


determining the population numbers of indistingulsliabic particles following Bosc-Einstein statistics 
using the grand canonical ensemble (Rubiucam, 1979). (Indistinguishable particles are chosen 
since the interchanging of particles does not affect tlie density distribution, which is the topic 
under discussion.) Tlie solution can thus be carried out in the usual statistical mechanics fasliion 
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(Morse, 1969, pp. 316*3I9‘, Reif, 1965, pp, 346*349), (1) is maximized subject to the constraints 
of the data (6) and S l*| * 1 : 


a_[- £ P, In Pj + a, P, p, F, , + S P, , ] 

dP, 


«o 


where , «q arc Lagrange multipliers, yielding Pj - exp(S where 




Z = c 0=pxp(£«,F,,) 


(7) 


Is the grand partition function. Using (5) in (7) gives 


Z = ^ exp [ctj S Ujj f| (tj ) + . . . + Oq Liijj fq (rj)3 


( 8 ) 


for Z, where ma^j has been redefined asaq. Note that the partial derivative of the logarithm of 
(8) with respect to fj ) gives 


ainZ 


d(«if,(rj)) 


f *’ji, 1 ? ”ji ^ 1 ^ *j ^ *^ • * ’ **■ *Q ^ *j ^3 


= S MjjPj “ '^Uj> 


(9) 


which is the expectation value of the number of particles in the cube at position Tj, This result 
will be used shortly. 

If there are no limits to the number of particles occupying each cube, then (8) can be factored 
as (Morse, 1969, p, 326; Reif, 1965, p. 347) 


z = nzi 

i J 
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where 


(fj)] iij} (10) 

tij«0 

M . V 

But this has the form x‘* , which is equal to 1 /(I -x). Hence (10) becomes 

n«0 

^ ^ 1 

J “ 1 -exp [o£j fj (^ ) -i . . . + «Q fq (rj)] (11) 

Therefore 


by (9). If it is now assumed that <iij> « 1, so that the Bose-Einstcin statistics pass over to 
Maxwcll-Boltzmann statistics (Reif, 1965, p. 352; Morse, 1969, p, 329), then the exponential 
term in tl»e denominator of (12) is much greater than 1, and (12) simplifies to 

<nj> « ‘ 

Multiplying this by m/dv and dropping the subscript j gives the expression 

„(■?) = ^ <.>= J2L + + (13) 

• dv J dv 

as the information theory density distribution. This is a most important result; it gives the form 
of the information theory density distribution for particles following Maxwell-Boltzmann statistics 
where the constraints on the density distribution have the form (2), The variance of the distribu- 
tion is discussed in Appendix A. 

The m/dv appearing in (13) is a troublesome factor. No commitment has been made either to 
the value of m or dv; and there seems to be no clear guidance on how to choose their values. As it 
turns out, this problem may be avoided by absorbing the factor into the spherically symmetric part 
of the density distribution, as shown next. 
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The Cask now is to simplify ( 1 3). Tliis involves two assumptions: first, that f| (r) pertains 
solely to the spherically symmetric part of the density distribution (expIaineU below); and second, 
that 

la2f2(ir) + ... + aQrQ(*r) 1«1 

so that 


0£»|f*,Cr) 't' * , ♦ 'i* i 

e-^ ^ ^ 1 +a2f2(r) + . . ,+aQfQ(r) 

where the fjrt, • . ► * fq (^) are tlie appropriate functions (4) for the gravity field (i.e, the 
anomalous density variation is small compared tlo the radial density.) With these two assumptions 
(13) becomes 


- ,-»'v m ®jfi(r) 

p,(r)a ^ c * * 

‘ dv 


I 4. ^ 

I *J* MM* I. ■■■ ■ « i 

(2fi+l)Mp4 


which may be written 


m (2£+l)Mpa® 


(14) 


where 




(15) 


with Pg(r) being the spherically symmetric part of p(r) and where subscripts find have been 
substituted for subscript rj. Obviously taking fj(7) to be 


f,(t) = — In [p„(r)dv/m] 
‘ a, ° 
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results ill (15). This is artlficiiil, to be sure; no value in (2) Is known for the earth where fj (t) 
has the form given above. However, it lias two advantages: the troublesome factor m/dv disap- 
pears, and any desired p^ir) can be used In (14), The Po(r) for the earth is known to a liigh degree 
of accuracy from other data: it cannot differ greatly from the PEM p^ir) of Dziewonski et al„ 
(1975). Hence it will be assumed here that the integral F, is known for the earth where fj(“r) is 
given by the above equation, in order to use Dzlcwonski ei al.'s (1975) nulial density distribution, 


All that remains to find the information tlieory density distribution is to evaluate the au,,^j 
using the gravity field coefficients. The raw gravity field coefficients given by (2) will not be used, 
however, for two reasons. The first reason is that the contribution of the earth’s hydrostatic 
equilibrium rotational bulge to the gravity field must be subtracted out, The second reason is that 
the gravity field of the ciwi umst also be subtracted from the coefficients when a crustal model 
is used, in this case the gravity field data become 


ft p;GEM -CR -HR 

~^fioi ”^fioi 


for the fi = even, m « 0, i = 1 terms and 

^ -GEM -CR 

for the other coefficients. The superscripts GEM, HE, and CR stand for "Goddard Earth Model,” 
"Hydrostatic Equilibrium,” and "Crust” resj)cctively. Actually only the £ = 2 and £ = 4 hydro- 
static equilibrium coefficients computed by Nakiboglu (1979) for Dziewonski et al.’s (1975) 

PEM will be used here; the higher degree hydrostatic equilibrium terms are assumed to be zero. 
Strictly, these terms should be included, but their computation Is difficult and the error in ignoring 
them is probably small. All of the crustal terms up to and including degree and order 36 will, 
however, be subtracted from the GEM lOB coefficients when a crustal model is lused, 
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Substituting (14) for pt"?) «nd for Cg,„j into (2) to evaluate the yields 

P,(R)»P„(R) [l+,S6se,„,R«Vj„„(9A)] (16) 

as the information theory density distribution, where the 


Sjg * 


(2g-H)PB 

3j^'^‘Jpo(R)R2«+2dR 


(17) 


arc found by using tSie orthogonality properties of the Vj„jj(0,X), where the earth is assumed to 
be a sphere, and where the variable r has been replaced by R * r/a^ for convenience so that 
0 < R < I . Also, pg Is the average density of the earth and Og « 637 1 km is the radius of the earth. 
Ry is the radius of the sphere in which the unknown density distribution to be inferred resides 
(the subscript u standing for “Upper.”) For example, Ry = 1 if no crust Is stripped off the earth 
and Ry - 6341/637 1 if a 30 km thick enist is stripped off, The integral in llie denonwnaior of 
(17) can be evaluated analytically, since Dziewonski et al. (1975) break up the earth into eight 
shells, with Pq(R) being givon as a polynomial in R in each shell. Table 1 gives the resulting 
from this computation for Rjj« 1 and tbr Ry » 6341/6371, 

RESULTS AND COMPARISONS WITH OTHER STUDIES 

The fundamental equations of tills paper arc ( 1 6) and ( 1 7). What they give is, in a sense, the 
broadest possible density anomalies: more localized anomalies arc not warranted by the data. 

Some general features of these equations should be noted before examining specific density distri- 
butions. 

The information theory density distribution p,(lr) (the subscript I standing for “Information 
Theory”) given by ( 1 6) is obviously a spherical harmonic expansion of the form 

P,(5) = P„(R) + Ap, =p„(R) + Epi„;(R) Yj„,,(e,X) (18) 




* 


10 








( 


i' 

k 
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ft 


where 

^P, ■ ^ ?in,|(R> Y«m|W.W. P«ml(R) " h R'^o<R>' ( 1! 

Note timt Ap, changes dlscontinuously when p^CR) does. Also, (16) Jias a form similar to that 
of tile equation giving a density variation due to lateral temperature differences (e.g, Phillii)s and 
Lambeck* 1980, p. 32)t 


P«Po n-«(T-To)J (20) 

Here a is tlie coefficient of thermal expansion, Pp and are the reference density and reference 
temperature, respectively, while T is the temperature, ([’he reason (16) has tliis form is due to 
the assumption that Ap| is small compared to pp( R)). So (16) is consistent with the idea that the 
density anomalies are due to lateral temperature differences- but does not necessarily imply that 
the anomalies are so caused. 

it should further be noted that there is no point in computing the power spectrum of the 
anomalous potential V| (AU) (e.g. Phillips and Lambeck, 1980, equation 1 1) generated by the 
density distiibution (16), since the Held coefficients are the given data, wliich (16) automatically 
satisfy. Moreover, the information theor>' density distribution (16) does not have a White noise 
spectrum, which Lambeck (1976) (see also Phillips and Lambeck, 1980) found will reproduce the 
observed anomalous potential spectrum (i.e. Kaula*sj[1967) rule-of-thumb). Instead the informa- 
tion theory density distribution Pj( R) concentrates the density anomalies towards the earth’s 
surface, due to the R® behavior of p^j^j(R) and the beltavior of 5^ in (19). (The 1!^ behavior 

may be seen by substituting Pg forpp(R) in (17) and evaluating the integral.) Also, it Is clear from 
(19) that the lower degree anomalies arc spread more evenly throughout the earth than the higher 
degree anomalies. This concentration of density anomalies towards the surface is in contrast to 
the findings of the Monte Carlo studies of Kaula (1977) and the mass-point studies of Lowrey (1978), 
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who indicate that the anomalies may increase with depth. Dziewonski et al. (1977) and Julian and 
Senguptu (1973), among others, also indicate that large anomalies are to be found deep within the 
mantle on the basis of seismic travel time studies. The statistical gravity study of Khan (1977), 
however, places most of the anomalies in the upper mantle, The seismic studies of Romanowicz 
(1979) and Cara (1979) and others show considerable upper mantle lateral structure, wlUle the 
analysis of satcllitc-to-satellite tracking data by Marsh et al. (1981) indicates that most of the 
gravity anomalies in the Pacific can be explained by lithospheric sources. 

It should be mentioned that the deiwity anomalies given by (16) and (17) are not confined 
to the crust and mantle, but extend into the core as well. To exclude them from the core the lower 
limit of the integral in (17) would have to be replaced by (the subscript L standing for “Lower”), 
where Rj^^ ~ 3485.7/6371,0, the radius of the core being 3485.7 km. In practice excluding the 
density anomalies from the cce makes little difference in the resulting density distribution for the 
crust and mantle, Finally, the boundaries where R) changes discontinuously are assumed to be 
spherical, so that there is no possibility of density anomalies arising from bumps on these bound- 
aries in the manner of Hide and Horai (1968) and McQueen and Stacey (1976), for example, 

The above constitute the general remarks on the information theory density distribution. 
Specific examples are discussed next, 

A computer program was written to produce maps in order to examine specific information 
theory density distributions.''' All of the maps are based on the GEM lOB gravity field (Lerch et al., 
1981). The GEM lOB field is based on satellite, surface gravity, and GEOS-3 altimetry data. Since 
GEM lOB is complete to degree and order 36, the limit of resolution is about 5 degrees of arc, or 
about 550 km on the earth’s surface. All of the higher degree terms (£>16) are assumed to be 
meaningful, although Pliillips and Lambeck (1980, p. 44) warn that these terms may largely be noise. 


*See Appendix B for a listing of the program. 
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Plate 1 shows the density variation iip, given by tl‘^) overlaid on (he global teetonie and 
volcanie aethlty map ol' bowman 1 1 ‘)8 1 ). The density variation is given on the surface of a sphere 
with radius 0368,0 km. It Is at this depth, 3 km, that (he oceans leave off and the rock surface 
begins in l)?iewonskl et al.’s Pl’M. No crust has been stripped off. The interval between 
contour lines is 0,002 g cm"^. This map looks quite similar to tlie CIEM lOB tVce*uir gravity anomaly 
map tS. Klosko, private communication, i980). U shows such typical featut'cs as the lows at 
Hudson Hay, Fennoscandia and some of the abyssal plain, s te.g., .Somali, Ilatteras); and higlus at 
some slow-moving ocean ridges te.g,, Mid-Atlantic, Sout Invest Indian Ocean), subduct Ion 70i>,es 
te,g„ Peni-Chile. Tonga-Kerma<loc), and hot spots te.g,, Hawaii, Iceland). Hence for qualitatively 
I'elating density and tectonies, similar to relating gravity to tectonics as done by Kaula (1972), a 
free-air gravity anomaly map mighi just as well be used. 


The map shows regions of artificial densities, due to the assumption that the earth is a sphere: 
the topography has been flattened, So at the Mkl-Atiunlic Ridge south of Iceland, for example, 
the effect of topography more than cancels the effect of the low density material upwelling beneath 
the ridge (assuming the basic validity of plate tectonics), producing a positive anomaly (Lambeck, 
1972) and hence an artillclally high density. 

Plate 2 shows the density distribution at 30 km depth (r "* o341 km), where 30 km of Iso- 
statically uncompensaCetl crust has been st ripped off the earth (giving essentially the Houguer 
anomahes). The interval bolsveen contour lines (s 0,02 g cm"^, Note that It gives low densities at 
some subduction xones. Removing the Andes, for example, completely erases the positive anomaly 
of Plate I, so that low densities prevaU. There is no sign of a highdensity subducting slab (wlilch 
is pmbably loo small to be seen in any case with the resolution employed here). 

Figures 1 . 2, and 3 show Apj on a plane which slices through the center of the earth In the 
equatorial plane for the cases of no ernst removed, 30 km ofisostatically compensated crust removed, 
and 30 km of isostatically unco«npensated cnist removed, respectively, The equatorial plane was 
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chosen because it illustrutes typical features of such slices, plus one atypical feature wlUch ap|«ai^ 

In Figure 2. 

All three figures illustiute a remark made earlier; that the Informatton theory density anom* 
alie^ cKtend deep into the earth, but the greatest density variation occurs near the surface. This 
is in qualitative agreement with Arkani-Hamed’s (1970) minimum shear strain energy density 
distribution, but his density anomalies are much lai-ger than those shown lieiw The figums also 
show regions where Apj changes sign with depth. This is also in qualitative agreement with 
Sanchez’s (1980) density distribution which minimizes the sum of the mantle shear strain energy 
plus gravllatlonal potei\tial energy. The size of the density variation shown in Figures I and 2 
is in fair agreenrent with Sanchez ( 1 980). However, the density distribution across Sanchez's (1980) 
slice through the equatorial plane looks notliing like those shown In the figures. It should be 
mcjitioned that the minimum energy solutions of Kaula (1963), Arkanl-Hamed (1970), and Sanchez 
(1980) give nonhydrostatic stresses which pi-obably exceed the finite strength of the mantle, Indi- 
cating that the assumed elastic rheology is unrealistic (tambeck, 1976, p. 6333), The alternation 
of tire sign of Ap, with depth is also in qualitative agreement with Lewis and Dorman ( 1970), who 
used a communications theory approach to the relation of density to topognqdry, 

The information theory density distributions .shown in Figures I to 3 clearly tend to form 
“pockets” of high and low density extending downwards from the surface. Tliere is no obvious 
convection pattern shown in any of the figures. None of them sirow what look like classical con- 
vection cells; that is, organized columns of low density material moving upwards and columns of 
high density material moving downwards. Pockets which slant at an angle to the local normal 
(such as the low density region at 305 degrees east longitude shown in Figure 2, for example) 
look like they might indicate some sort of liorizontai as well as vertical motion of material. But 
closer examination of the region around such features generally reveals that there are other pockets 
of similar density slanting towards tl\em (as is obvious with the two low density pockets between 60 
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and 90 degrees cast longitude shown in Figure 1, for example) which are not contoured. Hence 
the apparent “motion*' is probably an artifact of the contouring process, 

The one atypical feature mentioned earlier is the high density “blob” located at 285 degrees 
east longitude sliown in Figure 2, Closer examination of this featum shows that it is a “tube” of 
high density material connecting the Feru-Chilc Trench with the Middle America Trench. This 
feature is atypical in that the greatest density contnist occurs beneath the surface, while with the 
pockets the greatest density contrast occuis at the surface of the sphere inside which the density 
distribution is inferred, 

Figures I and 2 are quite similar to each other, The reason for this is that the removal of a 
30 km*thlck isostatically compensated crust affects mostly the high degree spherical harmonic 
terms, and not the low degree terms considered here. The density distribution shown in Figure 3 
is of course dominated by the removal of the uncompensated topography, and not the GEM I OB 
gravity field. 

If the density smomalies of Figures 1 and 2 are assumed to be due to lateral temperature 
variation, then typical temperature differences of 100 K are required, assuming cv in (20) is about 
3 X 10'^ K"^ (e,g„ Phillips and Lambeck, 1980, p, 53), The temperature differences must be 
about a factor of 10 greater to explain the density anomalies of Figure 3, assuming the same value 
fora. 

Maps showing the density distribution with the degree of the gravity and crustal fields restricted 
to fi < 1 6 also give pockets similar to those shown in the figures, Hence it appeare that the qualita- 
tive behavior of the information theory density distribution will not change if terms of higher degree 
(J^ > 36) than those considered here are included in the fields. 
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DISCUSSION 

One question which ariws is which of the three density distributions considered here Is to be 
prelerrcd, Obviously the one In whicli no crust is removed Is over-abstracted t no account is taken 
of the topC'fraphy. So it is not the preferred density distribution. The density distribution in which 
30 km of uncompensated crust is removed has intriguing consequences for the deep structure of 
continents (e.g., Jordan, 1975; but sec Anderson, 1979); very deep indeed, However, the geo- 
physical evidence favors the density distribution In which 30 km of compensated crust is removed; 
so it is preferred. But it is not ideal: it assumes that the topography is isostatically supported 
everywhere over the earth; and, moreover, with a depth of compensation of 30 km. This is certainly 
not the case. To cite just one e.vamplc where the assumptions fail, no account is taken of thermal 
isostacy at the ocean ridges (Ha.'cby and Turcottc, 1978), whcie the cft'cctive depth of compensation 
is greater than 30 km. Hence the density distribution still gives a high density for t!»e Mid-Atlantic 
Ridge south of Iceland, for e.Kamplc, after the crust is removed, when it should be low density. 

Hence better models of the crust are needed in order to use ITI to infer the density distribution 
below it, 

Another question which arises is why the information theory density distribution disagrees 
with the results of Lambcck (1976), Kaula (1977), and Lowrey (1978) which are also based solely 
on the relationship of density to gravity (and not on seismic travel times), and which give large 
density anomalies in the lower mantle. The answer appears to be tluit these studies examine only 
a limited number of models which mimic the external gravity field, while the information theory 
density distribution is a weighted average over all possible density models and reproduces the 
external gravity field exactly, 

The reason why the information theory density distribution disagrees with the scisnuc evidence 
of Julian and Sengupta (1973), Dziewonski et al. (1977) and others, which also give large anomalies 
in the lower mantle, seems clear enough: the seisinio data have not been included in the information 
theory approach. Their inclusion presumably would show large anomalies in the lower mantle. 
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The dominant impression from the foregoing remarks on the problem addressed licrc is one of 
simplicity. More data must be introduced to obtain results which are closer to the actual state of 
affairs inside the earth. Taking the earth to be made up of indistinguishable particles following 
Maxwell'Boltzmann statistics is an obviously simplifying, if fundamental assumption which must 
also be dealt with in order to obtain more realistic results. 

Introducing seismic travel times into the approaci) appears to be an obvious next step to take 
in examining the density anomalies. However, including the seismic travel time data into ITl 
appears to be difficult. It is not obvious how to put such information into an approach which sums 
over all possible models. Numerical models arc out of the question; dividing up titc earth into 10 
volume elements eacli of wliich may be occupied by up to 10 particles gives 10^^ models to 
consider — far too many already for a computer. Thus the problem must be done analytically, But 
even such a seemingly simple task as using free oscillation periods to obtain a spherically symmetric 
earth model has analytical difficulties; the data do not have the simple form (2), and the elastic 
parameters vary as well as the density, These difficulties make the evaluation of the partition func- 
tion Z troublesome (Graber, 1977). The variance (see Appendix A) and clioosing in/dv in (13) also 
post difficulties in applying ITI, 

On the positive side is the heart of the method; ITI (MEM) minimizes subjectivity. To illus- 
trate, there is no need in ITI to decide which lower degree harmonics to Ignore in examining density 
anomalies in the lithosphere (e.g., Marsh et al,, 1981), which is a highly subjective procedure. ITl 
weights all the harmonics automatically. And since ITl gives the one “best” (i.e., least subjective) 
model, it may well have use in choosing a Standard Earth Model, although the mathematical 
obstacles mentioned earlier would have to be overcome. 
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Plate I. Lateral density variation Ap, at 3 km depth (no crust removed). Note that no zero line is shown. 

which would make the map too busy. 





Plate 2. Latei^l density variat! >n Ap, at 30 km depth for 30 km of uncx>mpensated crust removed 
Note that no zero line is shown, which would make the map too busy. 
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FIGURE CAPTIONS 

Plate 1, Lateral density variation Apj at 3 km efepth (no crust removed), Note that no zero line is 
shown, which would make the map too busy. 

Plate 2. Lateral density variation Ap, at 30 km depth for 30 km of uncompensated crust removed. 
Note that no zero line is shown, which would make the map too busy. 

Figure 1 . Lateral density variation Ap| in the equatorial plane of the earth for the case of no crust 
removed. Sea level is displaced from the earth’s surface for clarity, and the rock-equivalent 
topography is greatly exaggerated, Anomalies in the core are not shown. 

Figure 2. Lateral density variation Apj in the equatorial plane of the earth for the case of 30 km 
of compensated crust removed. Sea level is displaced from the eartli’s surface for clarity, and the 
rock-equLvalent topography is greatly exaggerated. Anomalies in the core are not shown. 

Figure 3. Lateral density variation Ap, in the equatorial plane of the earth for the case of 30 km 
of uncompensated crust removed. Sea level is displaced from the earth’s surface for clarity, and the 
rock-equivalent topography is greatly exaggerated. Anomalies in the core are not shown. 


Table 1 


The Coefficients 6g for No Crust Removed (Ry » 1) and 
30 km of Crust Removed (Ry * (>341/6371), 


Degree fi 



c 

II 

Ry « 6341/6371 

1 

1,096 

1.112 

2 

2.794 

2,857 

3 

5,285 

5.445 

4 

8.580 

8,911 

5 

12,701 

13.295 

6 

17.659 

18.635 

7 

23.467 

24.965 

8 

30.130 

32.316 

9 

37,653 

40.715 

10 

46.036 

50.191 

11 

55,283 

60.771 

12 

65.392 

72.479 

13 

76.367 

85.344 

14 

88.207 

99.393 

15 

100.913 

114.653 

16 

114.486 

131.152 

17 

128,929 

148.920 

18 

144.241 

167.987 

19 

160.426 

188.382 

20 

177.485 

210.138 

21 

195.421 

233.286 

22 

214.235 

257.858 

23 

233.931 

283.890 

24 

254.510 

311.414 

25 

275.976 

340.466 

26 

298,332 

371.082 

27 

321.579 

403.298 

28 

345.722 

437.152 

29 

370.763 

472.681 

30 

396.706 

509,925 

31 

423.552 

548.922 

32 

451.306 

589.713 

33 

479,971 

632,339 

34 

509.548 

676.842 

35 

540.043 

723.263 

36 

571.458 

771.647 
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APPENDIX 


Finding the variance o,,^ » <nj^> - from (8) closely follows the staiulard statistical 
mechanics treatment (e.g., Reif, 1965, pp. 336-337), It is given by 



a^lnZ 


which gives 


<„j> 



m 

dv 


for particles following Bose-Einstcin statistics (Reif, 1965, p. 346; Morse, 1969, p, 333). Because 
it is assumed that (Uj) « 1, all that can be said about is that Op > p(j), since no commitment 
to the value of <nj> has been tuadc. Hence the data do not greatly constrain the density distribution. 


k 
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APPENDIX B 


Tht computer prokrtm ueed to generate the platea and (igurea ia Hated here. Phttea i and > uae 
aubroutlne VANDEIt. Flgurea 1, 2, and 3 uae aubroutine SLICE, 


•INFCKNATION rHECRV (ENSITY OISTAieUTION* 


PAI 


TPIS PRCCaaP COPPOTIS TtE INFOAPATION THtOAV 0EN$lt» OlStAIBUTtCN PAI 2 
EASED Oh IHE GRAVIfV FIILC. IT ASSURES NAXHELL-RCLTERAhN STATISTICS RAI ) 
AND THAT THE SPRCRICALU STANETAIC PART Of THE OENSITV OISTAIRUTION RAI 4 
IS KNCWh. pai ) 


THE PROCRAH USES SPHERICAL HARHOHICS. THE LATERAL OENSIIV 
VARIATtCh IS ASSURED TO HE A SMALL PERTURBATION ON TCP Of THE 
SPHERICALLY SYMMETRIC OINSIIY DISTRIBUTION. 


MAI a 

MAI T 

RAI a 


THE PROCRAM COMPUTES THE DENSITY VARIATICNt RELATIVE DEASITVt OR MAI 4 
ACTUAL DENSITY INSIDE THE EARTH. bHICH MAY RE DISPLAYED IN VARIOUS RA| 10 
MAYS ACCCROINS TO A CHOSEN SUBROUTINE. THE DISPLAY IS IN ONE Of TRO MAI U 
ROOESt PRINT-OUT PICTURES OR DICOREO PHOTOGRAPHS. RAI 12 

THE DENSITY IS OEfIREC AS THE SPHERICALLY SYMMETRIC MART RAI 1) 
IRHC SUB-EERDI PLUS THE DENSITY VARIATION lOELTA RHO). THE RELATIVE MAI I* 
DENSITY IS CEfINED AS DELTA RHO DIVIDED BY RHO SUB-ZERO. MAI 1$ 


CRUSTAL MODELS (SUCH AS AN ISOSTATICALLY COMPENSATED CRIST OR RAI 
AN ISOSTATICALLY UNCOMPENSATED CRUST I HAY BE SUBTRACTED Off If MAI 
DESIRED AND THE INFORRATiCN THEORY ALGORITHM APfLIEO TO THE REMAINING MAI 
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IT 
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GRAVITY SIGNAL. NAt 

ALSO. THE DEGREES USED BY THE SPHERICAL HARMONIC FIELDS CAN BE MAI 20 
RCSTRICIEC TO A CERTAIN RANGE IBEIMEEN LMIN AND LMAK INCLUSIVE* If MAI 21 
SC DESIRED. NAI 22 


NOTATION 


MAI 23 


GENERAL DATA MAI 2A 

PI • I.141592693SB9T900 MAI 2S 

PI0IV2 • Ml/2.000 NAI 26 

PI2 • 2.000*PI MAI 27 

LLOM • MINIMUM DECREE Of GRAVI TY/TOPOCRAPHY/CRUSTAL FIELD MAI 2B 

RECORDS READ MAI 29 

LUP • MARIMUM DECREE Of ORAVI TY/TOPOORAPHY/CRUSTAL FIELD MAI 30 

RECORDS RE/0 MAI 31 

LMIN • MINIMUM DECREE Of GRAVI TY/TOPOGRAPHY/CRUSTAL FIELD MAI 32 

TO BE USED IN PROGRAM MAI 33 

LMA* • MAXIMUM DECREE OF GRAVI TY/TOPOGRAPHY/CRUSTAL FIELD MAI 36 

TO BE USED IN FROGRAM NAI 3S 

SUBNME > NAME OF SUERCUTINE CALLED MAI 36 

COMENT « GENERAL COMMENTS ABOUT PARTICULAR RUN (SUCH AS THE NAME MAI 3T 

OF THE REGION BEING EXAMINED. ETCI MAI 31 

NMRINT ■ I IF GRAVI TY/TCP06RAPHY/CRUSTAL FIELD CCEFFICIENTS ARE NAI 39 

TO BE LISTED. *0 IF NOT MAI 60 

EARTH DATA MAI 61 

RHOAVC • AVERAGE CEMSITV OF THE EARTH IN CRAMS/CM**3 MAI 62 

R6DAVC ■ AVERAGE RACIUS OF EARTH IN KILOMETERS NAI 63 

RCORKM . RADIUS OF CUTER CORE IN KILOMETERS MAI 66 

RCORE « RCORKM/RAD/VC MAI 6S 

RLOMKM ■ LOMER RACIAL LIMIT ALLOhEO FOR OENSITV ANCMALIES. IN M6| 66 

KILOMETERS MAI 67 

RLOliER • RLONKM/RAOAVC HA| 61 

RUPKM < UPPER RADIAL LIMIT ALLOMED FOR DENSITY ANOMALIES. IN MAI 69 

KILOMETERS MAI 50 

RUPPER • RUPKM/RACAVC MAI 5| 

GRAVITY FIELD H6I 52 

GRAY • DESCRIPTION CF GRAVITY FIELD M6I 53 

CRUST Mat 56 

CCRUST « DESCRIPTION CF CRUSTAL MODEL MAI 55 

NCRUST • I IF CRUSTAL MCOEL FIELD IS SUBTRACTED FROM GRAVITY MAI 56 

FIELD. -C IF NCT MAI 57 

TCPCCRAPHV MAI 58 

TTOPC > DESCRIPTION CF TOPOGRAPHY MAI 59 

TOPC - ARRAY IN MUCH TOPOGRAPHY COEFFICIENTS ARE STORED MAI 60 

LATERAL DENSITY DISTRIBUTION MAI 61 

CONMNT « REMARKS ABCUT THE COEILI (SUCH AS THE DEPTHS FCR HHICH Hal 62 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


tHEY WWC tOfPirtD) t>*l 

CCE « INFOMXtTIOtflHEOftETIC CCEEFICiENrS EASED ON THE N*l 

SFHCAICAIL^ SYANETKIC DENSITY DISTlUflUnON NAt 

CLH ■ INFCANATIOA-IHiOAITIC AELATIVE OENSIIY/OENSIIY/OENSITY NA| 

VAAIATION IPFEFICAl HAAFONIC COEFFICIENTS NA| 

NCEN • 0 IF AELATIVE CENSITY DISTAIBOTION IS FECTTED. «l IF Ml 
ACIlAl density is FUTlEDf «2 IF DENSITY VAA|A||ON |$ Mi 
PLOTTED Ml 

HEFCEN ■ AEFEAENCE CCKSITYt USED ONLY IF NOEMl Ml 

PICTUAC DATA PA| 

NPHCTC • 0 IF PAINT'OLT PICTUAC IS OESlAEUt «t IF OICOAED NA| 

PHOIOOAAPII IS CE5IAEO PA| 

AAXSYA > NUA0EA OF ILPHANLHEAIC SYMBOLS LSEO FDA PAINT-CLT PAt 

PICTUAC CA OICCMED FHOTC MAI 

ALPFA « AAAAY IN UF|CH ALPHANUPEAiC SYMBOLS USED FCA PFINT»OUT PA| 

PICTUAES APE STOAEDf IN OAOEA FAOM LON TO HICH PAI 

AELATIVE OENSMY/DENSITY/OENSITV VAAIATION MAI 

BETA » AAAAY in WMCH DICOPED COLOA CO.OINC NUMBEAS AAE MAI 

STOAEDi USED FCA OtCQMEC FHOTO< IN OADEA FAOM LCW PA| 

TO HICH AELAIIVE DENSITV/PENSIIY/OENSITY VAAIATION MAI 

SCALE « AELATIVE OENSIIY/OENSITV/DENSITY VAAIATION INTEAVAL MA| 

BEThEEK ALFHANIMEAIC SYMBOLS lOA LDLOASI PAI 

lUNITSi NOPE IF NDEN<0, GKAMS/CP**S IF NOEN>lOP 21 PAI 


SUBRCUTINES CALLECt PAI 

1. CUTS (IF SUBNPE .10. CUTI PAI 

2. PIECE I IF SUCNPE .CD. FIE) PAI 

3. POLES I IF SUCNPE .EU. POL) MAI 

A. SECTNS I IF SUBNPE .EO. SEC) MAI 

5. SLICE IIF SUCNPE .E0> SLII PAI 

6. VANCER IIF SUBNPE .10. VAN) MAI 


NOTES PAI 

1. ALL SPHERICAL HARPOMC FIELDS ICRAVITV, TOPOGRAPHIC, CAUSTALI PAI 

USE KAULA'S A*P| NOFMALUATICN. PA] 

2. THE TOPQCRAPHIC FIELD IS N0APAU2E0 TO THE RADIUS OF THE PAI 

EARTH IIP KILCMETERE). PAI 

3. THE FOLLOWING CINCIIIONS PUST BE SATISFIEOl PAI 

LPIN .CE. LLOW PAI 

LPAP .LE. LUP . PAI 

A. THE DENSITY VARIATICN OR RELATIVE DENSITY IS MORE CCNVENICNT MAI 

TO CISPLAV THPN THE ACTUAL DENSITY. MAI 

S. IF WF HAVE ALFH«IJNt!Al WHERE JNUM .CT. MAXSVM. THEN THE MAI 

PAOCAAP! SETS JNIM^HAXSVM. ALSO, IF JNUM .LT. 1, THEN HE SET MAI 
JNUPi>l. LIKEHIST FOR BETAUNLMI. MAI 

IN OTHER HOADS, IF THE RELATIVE OENS ITY/OENSI TY/DENSI TY MAI 

VARIATION GOES EFF SCALE AT EITHER END, THEN THE LIMIT AT PAI 
THAT END IS USEC. MAI 

A. PROGRAM N0RPALI2ES ALL RADIAL DISTANCES BY DIVIDING RADIAL MAI 

CISTANCE UN KIIOPEIERST BY RADAVG. MAI 


SAMPLE INPUT 

DATA I COLUMN 1 

1 OF INPUT STARTS IN CCLUMN 3 HERE) 

MAI 

LLOH LLP 


IFORMATt AISI 

MAI 

36 36 



MAI 

LMIN LPAX 


IFORMATt AI5) 

MAI 

36 36 



MAI 

NCRUST NPRINT 


IFORMATt AI9) 

MAI 

1 1 



MAI 

RADAVG RCORKM 

RHCAVG 

IFORMATt BFIO.SI 

MAI 

6371.0 3A6S.7 

S.517 


MAI 



ICOPPNTI J) 

, 4*I,1A) 

IFORMATt 

t3A6,A2) 

MAI 

CENSITY 

COEFFICIENTS COEIL) 

(RANCEl 670 

TO 30 KM DEPTH) 

Mai 

rlowkm 

RUPKM 


IFORMATt 

BFIO.SI 

MAI 

970 UO 

63A1.0 




MAI 

K 

COEIK) 


IFORMATt 

I5,F20.Y) 

MAT 

1 

3.A909062 




MAI 

2 

6.AIC03B0 




MAI 

• • 

• 




MAI 

39 

723.66597 




MAI 

36 

T71.Y9AIY 




MAI 


IGRAVIJ), J*1,1A) IFORMATt ]3A6,A2I MAI 


A) 

AA 

AS 

AA 

AT 

Al 

AY 

70 

Tl 

T2 

TJ 

7A 

T5 

7A 

TT 

7B 

TY 

•0 

ai 

B2 

B3 

BA 

as 


BA 

BT 

BB 

BY 

YD 

Yl 

Y2 


YJ 

9A 

YS 

96 

YT 

YB 

YY 

100 

tot 

102 

103 

lOA 

109 

lOA 

107 

loa 

lOY 

no 


Ml 

112 

113 

llA 

MS 

MA 

117 

MB 

MY 

120 

121 

122 

123 

I2A 

129 
126 

127 

128 
12Y 

130 
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c 

CfHIOtf USING NlKlpDCLUlHTS) HVOliasunC FO« L«2iNt 

P>0 

PAI 

131 

u 

c 

INC S 

(FORPATI 

AX. 212, 2015. PI 

HAI 

132 

c 

1 0 0.0 0.0 



PAI 

133 

c 

1 t 0.0 0.0 



PAI 

134 

c 

• •• 



PAI 

135 

c 

}»}$ 0.]S«<OISTO-(I>O.TII97024I)-OR 



PAI 

13A 

c 

Itlt 0.2*«l«IIIC-ll>0.mM14TO>0« 



HAI 

137 

c 

c 

mopoiJi. 

(FORPATI 

13AA.A2I 

PAI 

I3f 

c 

DOCK'EOUtVHENr TOPCCNNrHY • NONMtUiEO TO 

RADIUS OF 

THE EARTH 

PAI 

159 

c 

c 

INC S 

(FORPATI 

AX, 212, 2015.1) 

PAI 

140 

c 

MCCET 1 0 C.T1IT20IIC-CN 0.0 



PAI 

141 

c 

MCOEE 1 1 G.0124S0R2C-CN O.N22EOS4JC>04 



PAI 

142 

c 

• • • 



PAI 

143 

c 

NECQEF)«3$ 0. l4S2iMS0>C9>0.7T«tSt ITC-ON 



PAI 

144 

c 

PECOEESt)* 0.72}SOai*C><»>0.24020MPC>0* 



PAI 

145 

c 

(CCPUSKJIi 

(FORPATI 

13AA.A2) 

PAI 

14A 

c 

CPPVITV COEEMCIENTS-UMCOEfENSATEO CRUST, FPON KAULPtLEE, SCR lEPS, NAVY 

PAI 

147 

c 

U N C $ 

(FORPATI 

2I4.1P2C19.A) 

PAI 

I4A 

c 

1 0 }.A392U0'CA 0.0 



PAI 

149 

c 

1 1 2.4S7P6SC-C4 2.072A7SD-CA 



PAI 

190 

c 

• • « 



PAI 

191 

c 

it if 3.S2UPP0-C7 «2.2AAA0I0'C7 



PA( 

152 

c 

3t 3t i.AA0A99D-Ca •I.U92AAD-C7 



PAI 

193 

c 

c 

NPHOTO NAASVN SCALE 

(FORPATI 

2I5.F10.5) 

PAI 

154 

c 

0 20 0.001 



PAI 

159 

c 

c 

(ALPHAIJI, J>l,2ll 

(FORPATI 

21A1) 

PAI 

ISA 

c 

JMiCFE0CPA12349A7P9f 



PAI 

197 

c 

c 

ISETAUl. J«l.llt 

(FORPATI 

11131 

PAI 

I9R 

c 

7 P 9 to 31 20 19 u 17 22 1 



PAI 

199 

c 

c 

NCEN PEFCEN 

1 FORPATI 

I5.F20.9) 

PAI 

lAO 

c 

0 0.0 



HAI 

lAl 

c 

• • • 



PAI 

IA2 

c 

IFOP the pest of THI SAFFLE INPUT, WHICH RELAICS TO THE 

PAI 

IA3 

c 

c 

PARTICULAR SU5ROU7INE CtUEC, SEE THE SUSSOUTINE.I 


PAI 

IA4 

c 

c 

IPPLICIT REAL*tlA>H,0>2l 



PAI 

IAS 


INTECER*2 RETAIIll 



PAI 

leA 


OIPENSICN ALPHAI23I ,LCUMI209.2I 



PAI 

167 


DINENSICN COPENT ll<l 



HAI 

lAR 


OIPENSICN CLPI3C,37,2I.C0EI36> 



PAI 

169 


OIPENSICN TOPOI36,3T,2I 



PAI 

170 


OIPENSICN COPPNTI IN l,CR AV( 14 1 ,CCRUSTI 14) ,TTQPO( 14) 

PAI 

171 


COPPCN/PLRP/CLP,CCE 



PAI 

172 


C0PPCNFeLPCFF,PI.PI!,PmV2,SCALE,NPHaT0,PAXSYP 


PAI 

173 


CCPPCN/PLKO/ALPHA.LCMHI 



PAI 

174 


COPPCN/PLRC/LP|N,LPAX,LLOH,LUP 



PAI 

175 


COPPCNFeLRH/RAOAVC.FCCRE,RLOPER,RUPP£R 



PAI 

176 


COPPCNFeLKlFCOPENT 



PAI 

177 


C0PPCN/eLKJ/NDEN,RE70EN 



PAI 

178 


CCPPCN/PLRP/TDPO 



PAI 

179 


COPPCN/PLKT/OETA 



PAI 

leo 


DATA POL.PIE/AHPCLE! .6FPIECE / 



PAI 

til 


DATA SLI.SEC.VAN/AHELICE .APSECTNS, AHVANDER/ 


PAI 

112 


DATA CUT/AHCUTS / 



PAI 

113 


DATA PLANK, DOLLAR/IF ,1FI/ 



PAI 

114 

c 

SUPPRESS UNDERFLOW ERROR PESSACES CN SUPROUTINE LCENOR 

PAI 

115 

(• 

CALL ERRSEri20P,2SA,>l, 11 



PAI 

186 

c 

initialue rasic constants 



PAI 

187 


PI«3.14I9924535P979E0 



PAI 

188 


F'PI/IPO.OOO 



PAI 

189 


PI2«2.0C0*PI 



PAI 

190 


PI0IV2*PI/2.000 



PAI 

191 


ALPHAI22I»00LLAR 



PAI 

192 

r 

ALPFA(23)>PLANK 



PAI 

193 

V 

c 

READ IN RASIC DATA 



PAI 

194 


REAC IS, 171 LLOW.LUF 



PAI 

195 


17 FCRPAT IAIS) 



PAI 

196 


REAC (5,171 LPIN.LPAK 



PAI 

197 


REAC 15,171 NCRUST.NPRINT 



PAI 

198 

L 

REAC (5,1041 RADAVC.RCOPKP.RHOAVC 



PAI 

199 


RCOPEoRCORKP/RACAVG 



PAI 

200 

c 

SET DATA arrays INITIALLY EQUAL TQ EERO 



PAI 

201 


35 



* t 


PAI 

202 


DC 21 L*l«Ll<M 


PA| 

201 


CCCIU*C*0U0 


PAI 

204 


OC 21 


PA) 

205 


icpnit»'tU»o.oco 


PAI 

206 


TCfCaiP»2)«0,QC0 


PA| 

20T 


CLPILiHtllsO.OUO 


PAI 

20 » 

ti 

l» 

ClP»liHi21»0.000 


PAI 

205 

c 

MAC IN CC^IU CMA 


PAl 

210 


«eAt 15,251 ICOPPNrUI, J«l,l4| 


PAI 

211 

ji 

fCHPAl U1A6,A2J 


PA| 

212 


At AC 15,1041 AUWAN.HOI'AH 


PAI 

211 

IfiAl 

fCAPM ISR0.5I 


PAI 

214 


HlOktA*AlOMKH/AADAV( 


PAI 

m 


AliPPtR>MUPKP/AACAVO 


PAI 

216 

ICCIM'O 


PAI 

217 


CC 24 L«LMN,t.PAA 


PAI 

2I» 


REACtS, 2U K.COEIRI 


PAI 

215 


FCRPM U5,f20.<JI 


PAI 

220 


ICCCM»ICquNI « 1 


PAI 

221 

24 

t 

c 

CCNIINUC 


PAI 

222 

RIAC IK CRAVliy FlHO 0‘S ANO 5«S 


PAI 

221 

c 

RfAC 15.251 IGRAVIJI, J<1.14) 


PAI 

224 


URlIt U,|4) 


PAI 

225 


HMU U,12) 


PA) 

226 

12 

fCRPAT llOK.'CKAVin fiao C0EfFltl£NtS*,/I 

PAI 

22T 


WRITE U.Jll tCRAVii.1, J«1,14I 


PAI 

22R 

11 

FCRPAT l20K,tlAA,A2.///| 


PAI 

225 


WRITE U,14l 


PAI 

210 

14 

fCRRAT IIOR,' L',2K,» H'.RH, <C* , 10*. 

•S', /I 

PAI 

211 


JCOL’KT*C 


PAI 

212 

2 

REAC t5il<£WC>}l L.K.L.l 


PAI 

213 

1 

fCRRAT IKX, 212. 2015. «l 


PAI 

214 


IF IMPRINT .EQ. 01 CU TC 111 


PAI 

215 

C 

WRI1E CRAVITV HELD C *S AKO 5'5 


PAI 

216 


WRITE Ifc.lSI t.P.C.S 


PA) 

217 

IS 

FCRPAT IIOX.U.2A.I2,2X.ni5.a.2X,Ql5 

.61 

PAI 

21« 

111 

JC»%KT=JCOUfiT t S 


PAI 

215 


IF IL .CT. IPAXl CO TC 11 


PAI 

240 


If IL .LT. LPINl CO TC 11 


PAI 

241 

c 

CCRPUTt DCNSITV OISIKieu ION COEFFICIENTS 

PAI 

242 


CCtrF»HfOAVC»ICCEIUI 


PAI 

243 


Rl» P ♦ 1 


PAI 

244 


ClWa,NI,lloC*CCEFF 


PAI 

245 


aPIL,Rl,21i5»COErF 


PA) 

246 


IF a .EO. LCP .AKU, P .EO. LUM CO 

TO 1 

PAI 

247 

11 

CCNTINUC 


PAI 

24R 


OC TC 2 


PAI 

245 

1 

CCNTINUE 


PAI 

250 


WRITE ItalO) JCDUI/T 


PAI 

251 

120 

FCRPAT l/.lOX. *NUP'ltR OF CRAVITV FIECO RECORDS READI • . IX, 1 101 

PAI 

252 

b 

c 

REAC IK TCPOGRAPHIC C*S AKD S*S 


PAI 

251 

U 

REAC 15.25) lITCPCai, 2*I.14I 


PAI 

254 


WRITE 16.141 


PAI 

255 


WRITE <6.1211 


PAI 

256 

121 

FCRPAT l20X,*TOP0CR»PhV COEFFICIENTS 

'./I 

PAI 

257 


WRITE <6,111 <ITQPO<JI. J>1,)4I 


PAI 

256 


WRITE (6.14) 


PAI 

255 


L€OCM»0 


PAI 

260 

122 

REAC t5,l,ENC024) L.P,C,S 


PAI 

261 


If IKPRINT .EO. 01 CP IC 112 


PAI 

262 

C 

WRITE TCPCCRAPFlC C'S AKO S'S 


PAI 

261 


WRITE 16,151 C.P.CtS 


PAI 

264 

132 

LCOCKT-ICOUNT • 1 


PAI 

265 

IF a .CT. LPAXI CO TC 221 


PAI 

266 


IF a .IT. LPINl CU TC 321 


PAI 

267 


PI*P • 1 


PAI 

266 


Tcpca,Pi,u«c 


PAI 

265 


TCPCa,Pt,2l-S 


PAI 

270 


IF a .EQ. CUP .XKU. P .EQ. >p) GO 

TO 324 

PAI 

2Tl 

121 

CCNTINUE 


PAI 

272 


GO TC 122 


PAI 

271 

124 

CCNT INUC 


PAI 

274 


WRITE <6,1251 UOliNT 


PAI 

275 

125 

A 

FCRPAT </,lOX,‘KUP8IR Of T0PO6RAPHV 

RECORDS R£AD<'rlX,l 101 

PAI 

276 

C 

IF IKCRUSr .EQ. 01 PCCUKT>0 


PAI 

277 

c 

IF IKCRUST .EQ. 01 CO TC 41 


PAI 

276 

c 

REAC IK CRUSTAL C'S AKO S'S 


PAI 

275 


c 
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» 


* 


HtAC iS«25) ICCKUSrUI* 

WMtE lEiU) Ml 

H«ITE ie.Z6) ICCALSIlJlf Mi 

KCOlNr»C Ml 

WRITE U.J4I WAI 

At REAC l9tAZ»fK0M3l 1|W>C«S WAI 

AZ ECAWAI IZIA,lf>20i;.(i WA| 

IF IKPRINT »E0, 01 (0 TC 333 WAI 

C WRITE CUT CAUSIAL C'S AK0 S«S WAI 

WAITE Iti33l LtW.Cil WAI 

333 ACDLM«KCOUNI • | WAI 

IF a .CT. IWAXI CO TC J2 WAI 

IF a .IT. LWlfO CO TC JZ WAI 

CCEFF^AFOAVCAICCEIL)) WAI 

Wt>W » I WAI 

c sunrRAcr crustal fieu each oeiisitv coefficients wai 

CLwaiHi,u«CLHa.Mi.n > occeff wai 

CLWa.Mt.2»>CLHa.H|.2l - $*C0EFF WAI 

IF a .EO. LLP .AND. W .EO. LLP) CO TO AT WAI 

SZ CCNTINUE WAI 

GC 1C At WAI 

A3 CCNTINUE WAI 

WAITE UtlZl) KC0LN1 WAI 

^IZl FCRPAT I/. lOA.aUPflfR OF CRUSTAL FIELD AECOAOS READl'.lAf 110) WA| 

C REAC IK display DATA WAI 

REAC IS.ZOtl NPHOTD.WAXSYW. SCALE WAI 

201 FCRWAT I2I5.TF10.S) WA| 

REAC IS. $01 lALPHAlJ). J>1.2U WAI 

50 FCRWAT (21A|) WA| 

REAC IS. SI) ieETAiJ),J>t.lt) WAI 

51 FCRWATnil)) WAI 

REAC 15.21) KDEN.RETDEN WAI 

i 

C WRITE CUT INPUT DATA WAI 

C 

WRITE lE.lAI WAI 

WRITE IE.300I WAI 

3C0 FCRWAT (20X.MNFaRM/TI0K THEORY DENSITY 01 STRIBUTICN* . //.20X . WAI 
. 'BASED ON THE GRAVITY MELD COEFFICIENTS' ./.20K. WAI 

. *AKC KNOWN SPHERICALLY SYWWETRIC DENSITY DISTRIBUTION' ,////) WAI 

WRITE IE.3A0) LLOW WAI 

340 FCRWAT IIOX. 'LLCW^ > .1 4. TOR. < I )«{NiWUW DEGREE OF SPHERICAl HARWaNIC WA| 

.r-IELL'SI'./) WAI 

WHITE It, 141) LUP PAI 

341 FCRWAT UOX.'LUP >' .1 4. 30K. ' I WAXIWUM DECREE OF SPHERICAL HARWONIC WAI 

.FtELCSI'.A/l WAI 

WRITE It. 3421 LWIN WAI 

342 FCRWAT IIOX, 'LWlN>'. 14. iOX.'IWINIWUW DEGREE USED IN RUN)'./) WAI 

WRITE ie>343! LWAX WAI 

343 FCRWAT I lOX. 'LWAX*' . 1 4. 30X. • I WAXIWUM DECREE USED INRUNI'.//) HA| 

WRITE It. 3031 NCRLS1 WAI 

303 FCRWAT UOX«'NCRUST''iI3.23X.«(»l IF CRUSTAL WOOEL IS SLBTRACTEDI WAI 

.«0 IF NCT)'.//) WAI 

WRITE lt.3301 NPRIN1 WAI 

330 FCRWAT IIOX. 'NWR|NT«' . I 3.29X. • (■ 1 IF GRAVITY. TOfOCRAPHIC. AND CRUWAI 
.STAl COEFFICIENTS AWE PRINTEOI *0 IF NOT)'.//) WAI 

WRITE <6.3041 RAOAVC WAI 

304 FCRWAT IIOX. 'RAOAVC*'. IX. FlO.S.lX.'KW. tax, 'lAVERACE RACIUS OF EARWAI 

.THI'.//I WAI 

WRITE It. 3051 RCORKW WAI 

305 FCRWAT I lOX. 'RCCRKW*' . I X.FI0.3. IX. 'KW' . IflX, ' IRAOlUS OF CDREI'./VI WAI 

WRITE lt»306l RHOAVC WAI 

306 FCRWAT IIOX«'RHCAVC«'.lX«Fia.5.1X«'CRAW/CN**3'.10X,'|AVERACE OENS|WA| 

•tY CF EARTHI •,//.) WAI 

WRITE It. 307) RLOWKW WAI 

307 FCRWAT COX. 'RLCWKW*' . 1 X.F10.3, IX, 'KW . IBX. • ILOWER BOUNC ON ANOWALWAI 

•|ES)'.//I WAI 

WRITE It. JOB) RUPKW WAI 

3CB FCRWAT ( lOX. 'RUPKW '.2X .F 10.3. IX .'KW . 18X. • lUPWER BOUND ON ANOWALIWAI 
.ES)'.//) WAI 

WRITE (t.309) NPHCTC WAI 

309 FCRWAT I lOX. 'NPHOTO*' . I 3. 29X. ' I * 1 IF PHOTO IS TAKEN! >0 FOR PRINT-WAI 

.OUT PICTURE!'.//! WAI 

WRITE lt.311) SCALE WAI 

111 FCRWAT IIOX. 'SCALE*'. IX. Ft0.5.22X.'IUNITSs NONE IF NOEN*0. CRAW/CWWAI 
.**3 IF NDEN*I OR 2)'.// I WAI 

WRITE lt.310) WAXSYW WRI 

310 FORWAT <IOX.'MAXSYW>'il5.2TX.'INLWBER OF ALPHANUWERIC SYWBOLS USEDWAI 

. )•.//! WAI 

WRITE It»3l2) (ALPHAIJ). J*l,21) WAI 

312 FCRWAT I lOX. 'ALPHAI Jl :• .IX. 21A1, BX.'IFROW LOW TO HIGH CENSI TY ) • ./WAI 

./) WAI 

WRITE It. 3131 IBETAIJI. 3*1.111 WAI 

313 FCRWAT I lOX, '6ETAI3 X '. IX, I II 13. tXI . 10X« • I USED FCR WHOTC)',//! WAI 

WRIIE It, 314) NCEN WAI 

314 FCRWAT ItOX.'NOEN*' ,I3.31X.'I*0 FOR RELATIVE DENSITY. *1 FOR ACTUAWAI 

.L DENSITY. *2 FOR DENSITY VARIATION)'.//) WAI 


210 

2Bl 

282 

Z«l 

284 

285 

286 
287 
218 
219 

290 

291 

292 

293 

294 

295 

296 

297 

298 

299 

300 

301 

302 

303 


304 

305 

306 

307 

308 

309 

310 

311 

112 

113 

314 

315 
116 

317 

318 

319 
120 
121 

322 

323 

324 

325 

326 

327 

328 

329 

330 
311 

332 

333 

334 

335 

336 

337 

338 
139 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

152 

353 

354 

355 

356 

357 

351 

359 

360 
161 

362 

363 

364 


37 
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1 

-i 


Jlf 


Umif REFCC^ l>k 

Ep«»-*t OOII*'K{EOIN«SI>.EIO,7il)l,«SK*N/CH«*|>,tOl«M>lrEEE*tCI OfNM 



.$mi USED KBEN NCC*>1I'*//) 

PAI 

1*7 


MRITE |«,l7*t 

PAI 

!*• 

12* 

ECMAA7 «//«IOAl 

RAi 

1*4 


H*ITf «*«1I9) (CRAViJIi J>l»l*l 

RAI 

170 

111 

.FORMT UOKi'CRAVin MUDI *t |Xf 11A*>A2I 

RAI 

17 1 


NRIIE t*il20) JC0UN1 

RAt 

172 


URIIE U.12$l mCRCUl. J«M*» 

RAI 

171 

12$ 

FORRAT <//•l0X»•TCRC0RAFHY•«]X••iMX>llA6«A^) 

RAI 

17* 


RRlTf U.12$) LC0UN1 

RAI 

175 


IF IKCRUST .EO. 01 CO TC 31* 

RAI 

17* 


NRIff |»»ll*l ICCRUHIJI* J>UI4» 

RAI 

177 

lU 

FCRRAT l//il0X»*CRUiTS(X«M<ilXtl}**iA2> 

RAI 

17B 


URI1E Utlll) KCDUNI 

RAI 

179 


CC 1C 114 

RAI 

no 

11* 

MRITf UillOl 

RAI 

1*1 

120 

FCRRAT l//•lOX••CRUiT•.eX••l•tlXt•NO CRUSTAL RODEL SUBTRACTED OUT'RAI 

1B2 


♦ \ 

RAI 

111 

11* 

CCNf INUI 

RAI 

11* 

WMie CUT COfilU 

RAI 

115 


MRITE l*«l*l 

RAI 

11* 

I* 

FCRRAT UHM 

RAI 

1*7 


WRITE U.2*l ICORRNIUIt J<l«l*l 

RAI 

111 

21 

FCRRAt U0X,11A*||A2,//) 

RAI 

1*4 


WRITE l*»29) 

RAI 

140 

2* 

FCRRAT <2*X,<L*i2CX,<COE(LI*,/| 

RAI 

141 


DC 3C L>LM|N«IMX 

RAI 

142 

1C 

WRITE lA.in L.comi 

RAI 

141 

11 

FORRCr l20X»IS»IOXtF2C«9l 

RAI 

14* 


WRITE I*il22l ICOUNI 

RAI 

145 

122 

FCRRAT l//(9X«*KUFBtR OT COEILI RECORDS READt > , IXi 1101 

RAI 

14* 

RC*C IN SUBNOUTINE NAIE 

RAI 

147 


REAC i9t*9l SUBNRE 

RAI 

141 

*9 

FCRRAT |A6» 

RAI 

144 

REAC IK SbBKQUTINE C01KEK1S 

RAI 

*00 


REAC |S«2$) ICQREKTlJit J*l«l*l 

RAI 

*01 

DECIDE UHICH SUBHCUIIIE 1$ CALLED 

RAI 

*02 


IF ISUBKME .EO. SLI 1 CO TO IOC 

RAI 

*01 


IF ISUBKRE .EO. SECI CO TO 101 

RAI 

*0* 


IF ISUBKHE .EOf VANI CO TO 102 

RAI 

*05 


IF ISUBftME .Efi. FCLI CO TO 101 

RAI 

*0* 


IF I5UBME .EO. F|E» CO TQ 10* 

RAI 

*07 


IF ISURKRE .EO. CUT I CO TO 107 

RAI 

*0* 


CC TC 94 

RAI 

*09 

ICQ 

CONTINUE 

RAI 

*10 


REAC ISfAOOl PHIl,XlNCAliFH|2iXLROA2,TaLERD»NRAO 

MAI 

*11 

ACO 

FCRRAT (SFlO.SilSI 

RAI 

*12 


CALL SLICE !FHIl,XLFOAl.PHI2,*LROA2,TOLERO#NRAOI 

RAI 

*11 


GC TC 99 

RAI 

*1* 

101 

CCNTINUE 

RAI 

*15 


REAC 15.171 NSEC.KRAO 

RAI 

*1* 


CALL SECTNSiNSECiMAOl 

RAI 

♦ 17 


GC TC 99 

RAI 

*1* 

102 

CCNTINUE 

RAI 

*19 


REAC 15.2021 NRAO.RXR.HXAF.WNAP 

RAI 

*20 


CALL VANOER IRKf.HRAF.WRAF.NRAOl 

RAI 

*21 


GO TC 49 

RAI 

*22 

ICl 

CONTINUE 

RAI 

*21 


REAC 15.2021 NNAD.A*C1>A>,RNR 

RAI 

*2* 

202 

FCRRAT II5.1FI0.5I 

RAI 

*25 


NFOLEol 

RAI 

*2* 


CALL FOLES IRKR.ANCFAX.NRAD.NFOLEI 

RAI 

*27 


NFCLE«2 

RAI 

*21 


CALL ROLES (RXRiANCRAX.NRAD.NROLEl 

RAI 

*29 


CC TC 94 

RAl 

♦30 

10* 

CCNTINUE 

RAI 

*11 


REAC I5.**l RI.R2.rFETA|.IHCTA2.XLHDAl,XLNOA2.NSfACE.NSLRf 

RAI 

*12 

*A 
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*11 
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. NCLT.NNORE.XPACNI) 
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STOP 
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END 
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FINC IFE PEIAIIVE DtN:iiy« OENSirVt OR OEMSItV V«RIATI(W CN EACH 
HWIICM41 UNI 

CCNEUfE TFE RCUnvE CfNSlfV* OENSifyp OR OENSITV VRRIRriCN U ROIN 
rFt$ SUIROUTINE CCHFUIE! THE RELATIVE DENSIIVf 0EN$ITVt CR 
CENSnV VRRUriON CN VERTICAL RECIRNCLCS ICUTSI NHICH ARE 
RERPENOICLLRR TO A CREAI CIRCLE CCNNECTINC TW POINTS CN THE EARTH'S 
SURFACE AT INIERVALS EQIALIT SRACEC ALONS THE GREAT CIRCLE lECRENT. 
FENCE NFEN LL'OMNC CORN ON tHE SURFACE OF THE EARTH IT LOCKS LIKE 
IHE LACINCS OF A FOCT0AU |F|IH THE CUTS ERTENOINO VERIICALLV 
CChNNARCSIi 


D C 0 0 

0 C 0 0 

OliMliXSXCRSRRXRORRRiiXRa 
0 C 0 0 

0 C 0 0 

CUTS I C 
GREAT CIRCLE I X 


EACF CUT IS CISFLAVIO ALONG WITH THE TOROCRARHV ALCNG TFE TRACK 
CF THE cur. 

ALL CUTS HAVE THE SANE LENGTH ANC DEPTH. ALL POINTS ON A VERTICAL 
PRfNT-OUT (PHOIOI LINE FAVE THE SANE LATirool AND LONGITUDE. 

TFE CUTS ARE OISPLATEC AS RECTANGLES* $0 THAT THERE IS A CERTAIN 
ANQUNT CF DISTORTION. HI ARE LOOKING AT THE CUTS IN IHE FCLIONING 
kAVt START NHERE THE CRTAT CIRCLE INIERSECIS THE ESUATCR ANC CO NORTH 
along TFE GREAT CIRCLE LNHI THE CUTS ARE REACHED. 

The INO POINTS COKNECIEC BV the great CIRCLE SeOPENT NAT BE 
CHOSEN AT WILL. THE NLNFER CF CUTS* THEIR LENGTH. AND THEIR DEPTH NAT 
AL5C BE CFOSEN AT WILL. TFERE IS ALNAVS A CUT AT EACH ENDPOINT OF 
IHE GREAT CIRCLE SECNENI. 


NCTATION 
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XLNCAI 

PFIJ 

XLNCA? 

OEPIFl 

OEPIFB 

ARCING 

rCLERC 


NCUT 

NFORE 

XNACNI 


LATIIUCE OF FIRST GREAT CIRCLE ENOPDIHr IN BteRfES 
LONGITUDE CF FIRST GREAT CIRCLE ENDPOINT IN DECREES 
LATITUCE OF SECOND GREAT CIRCLE ENDPOINT IN DECREES 
LONGITUDE CF SECOND GREAT CIRCLE ENDPOINT IN DECREES 
DEPTH TO LCWER CUT BCUNCARV IN KILONETERS 
DEPTH TC UFPER CUT BCUNCARV IN KILONETERS 
LENGTH OF CUT IN DEGREES 

TOLERANCE UNIT IN DECREES. BELOW WHICH A SNALl 
DIFFERENCE IN LAIITUCE CR LONGITUDE OF THE GREAT CIRCLE 
ENDFOINfS IS SET ECUAL TO EERO TO AVOID SINGULARITIES 
NUPIER OF CUTS ALONG THE GREAT CIRCLE SECNENT 
LENCIH OF cut IN SPACES ACROSS PACE lOR PHCTOI 
NACNIFICATION CF DEPTH CONPAREO TO LENGTH CP CLT lol.O 
FOR TRUE OIPTH TO LENCIH RATIOI 


SUBRCUTINES CALLECi 

1. ANCLE 

2. CENSrV 

3. TCPCCR 
A. VCC 


NCTES 

1. XLNCAI .LC. XLNCA2. 

2. CEPTHI .CT. DEP1H2. 

3. NhORE IS ODD ANC .LE. 121. 

A. CEPTHI IS USUALLY CFANCEO SLICNTLV BV THE SUBROUTINE TO HAKE 
TFE VERTICAL LENGTH OF THE CUTS AN INTEGER NUNBER OF SPACES 
long on the PRINT-OLT PICTURE (OR OICONEO PHCTCCRAPFI. 
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CUT 

29 

CUT 

2A 

CUT 

27 

CUT 

21 

cur 

29 

CUT 

30 

CUT 

31 

cur 

32 

CUT 

33 

CUT 

3A 
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99 
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99 
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cur 67 
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|66LICir 6E6L*«(6-H,0-lt 
l6rtCER«2 T6E6M2I). EEUmi 
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OIT'EKSICN K66PI121I 
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CCPPCN/eiXD/ElPM.lCUPJ 
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CCHPCK/EUI/COPENf 
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CCPPCK/EEXT/RETX 
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cur 
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cur 

cur 

cur 
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cur 
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cur 

CUT 

CUT 

cur 
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77 

71 
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cur S7 
cur IE 
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C FINC PICPIE or STPBCl/CCUCR RXNCE 
SVPMXxPAXSyP 

SVPIU«SVPMAX/2.0CC « C.500 

c 

R1*I.0DC - (OEPTHl/PACAVC) 

R2«I.OOC - (CEPTHE/PACAVCI 
HCM 16 HCAI - I 
XMCCL«(U0R772.00CI » l.OOO 
PICCUI-INHOA2 - 11/J 
$rHECH«5.0DO/16,OCO*XPA<Nn 
ir IKPKCTO .EO. II !rPECH>1.0CO/XHACNI 
rACrCA»HOA2/)R2*APClNa»}TPeCH*rl 
NAy>IP2 - RIIAPACtaP • 1.000 
Py»MY CUT 100 

Ri*x 2 - Ry/rAcro« cur loi 

DEPT»<1«PADAV0*I1.C0C - 111 CUT 102 

AAe6P'»ARcU»S*SADAVe*R2*f CUT 101 

C PRINT INPUT DATA CUT 106 

UPirC (6,711 CUT 105 

71 rCRPAT IIHII CUT 106 

WAITE 16,1001 CUT 107 

ICO rCAPAT l7//,10X,*SUfACU1INE CAUEDI CUT$«,///| CUT 101 

WAITE 16,1111 ICOPEMUI, J»I,I6| CUT 109 

111 FCAPAT I10X,11A6,A2,///I CUT IIO 

WAITE 16,671 PKIt,XlWCAl CUT 111 

67 FCAPAT (lX,'FiAST CAEAT CIACLE POINT I ' ,5X, 'lATITUDE-* , IX.FIO.S, CUT 112 

. 1K,«CECAEESMOX,UONCITUOE«',IX,F10.5,IX,«DECAEES',7/I CUT 111 

WAITE 16,661 PKI2,X1PICA< CUT 116 

66 FCAPAT llXt'SECCNC CAEAT CIACLE POINT t • ,5X, HATI TUOE>* , lX.FIO.S, CUT 115 

. IX,«CECAEE$*,IOX,'IONCITUOE>«,1X,FI0.5,IX,'DEGAEES',/7I CUT 116 

WAITE 16,651 0EPm,DEPIH2 CUT 117 

65 FCAPAT ISX.'DEPTHS FELOP SURFACE AETUEENMX,FI0,5,1X,*AP AND*, IX, CUT III 

. F1C.5,IX,*XP*,//I CUT 119 

WAITE 16,611 NCUT,NF0P{,XNACNI,AACLNC,AACAN CUT 120 

61 FCAPAT I5X,*NCUT>*, I5,1CX.*NHOA2>*,IS,10X,*XNACNI«*,1X,F|0.5»10X, CUT 121 

. •AACLNC-*,IX,FI0.5,IX,*DECA|E5 * • ,F 10.5, IX, *AN* , //} CUT 122 

NCEC>1 CUT 121 

RAOPAP«;6.0DO CUT 126 

WAITE 16,1121 AAOPAF CUT 125 

312 FCAPAT (7//, IX,*X,y CAID COOAOINATES ON VAN OEA CAINTEA HAP WITH CUT 126 

.AADIUS*.IX,FIO.S,IX,*CH*,//l CUT 127 

CALL VOC |PHI1,XLP0AI,M0HAP,N0EC,X,VI CUT 121 

WAITE t6,}lll X,y CUT 129 

113 FCAPAT (IOX,*FIAST CHEAT CIACLE POINT t • , IX, *XCA10« • , IX, FI. 1, IX, CUT 1)0 

.•CP*,10X,*yCRIDa*,lX,FI.),lX,*CP*,7/l CUT 111 

CALL VOC (PHI2,XLPO«2,MDHAP,NDEC,X,yi CUT 1)2 

WAtTE 16,3161 X,y CUT 1)) 

316 FCAPAT (lOX, 'SECOND CAEAT CIACLE POINTt * ,2X,*XCRI0>* ,1K,F«.),1X, CUT 136 

. •CP*,10X,‘VCAIC>',IX,FE.).1X,*CP*,//I CUT 1)5 

C CCNVEAT ANCLES FAOH OfCAEfS TO AACIANS CUT 1)6 

PF|l«PHll»F CUT 1)7 

PM2«PMI2»F CUT III 

XLHCAl«XLHDAI«r CUT 1)9 

XLHCA2>XLHDA2*F CUT 160 

AACLNCAACLNC'F cut 161 

rCLEA>TCLEAC*F CUT 162 

C FINC Ml, CHECA, COSXM, EINXXI, CCSOHC, ElN!;*<C CUT 16) 

TCLPKI«CSCAT(PH1I**2 a FHI2M2I CUT 166 

C IS THE CHEAT CIACLE TFE ECUATOAT CUT 165 

IF ITCLPHI - rCLEAl 79,10,10 CUT 166 

79 XXI'O.OCO CUT 167 
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32 

31 


55 

55 
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57 
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«5 
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59 
62 

61 
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c 

c 
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115 
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165 

StN»M«C.OOO 
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169 

01 >tC 6 > 0<000 
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CCSC 5 C-I .000 
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siNcn«o.ooo 
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PSmi«KLM 061 
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5 SII 2 t>XLHOI 2 


CUT 
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TCU 65 >C 65 $mM 062 - 3 L 7 D 61 I 


CUT 
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IS IHI 05167 CI 5 CLE 6 5 f 5 lO| 6 N 7 


CUT 

157 

If ITCU 6 M - TDLE 5 I I 2 tl 3 t 61 
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CCSXX(> 0 .ORO 


CUI 
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05 EC 6 -XLMC 6 I 


COT 

162 
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fS|( 2 W 5 HI 2 
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0 PEC 6 - 06 TANU 6 NC 5 Ct 
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CUT 
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If (TANXXII 52 f 92 i 53 


COT 
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COT 
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t 6 NMI» 7 ANXXI 


COT 
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XX|>C 6 T 6 N|T 6 NXXI| 
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CUT 
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CUT 
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cut 
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cut 
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cut 
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cur 

306 

6 

66API6UAI.PHAI19I 

CUT 

110 


WRI1E 11,161 

CUT 

1U 


HM1E 16,2061 

cut 

112 

t«6 

PC66A1 |6X,166S6»6,M0C1S INIIICAIE SEA EEVEUM1 

CUT 

111 


DC SCI LP*l,61 

CUI 

116 


U»6S “ CP 

CUT 

HI 


PU»Ct 

CUT 

H* 


lltlCH1«irU *. SI. tiKI 41.000 

CUT 

HI 


00 SCS 6«1,hPQN< 

cut 

HP 


a 


ir ut ;ii x»ahki*lhu 



COT 

IlY 

If UXUKtKl .BQ. lit XMI>lKI*tKS 



COT 

>20 

}QI ICMINUB 



COT 

>21 

wMtc HiictiiiilixptM) XKitian 



tui 

>22 

ici rcHXM iix»r4.i»iimi 



CU1 

>2> 

If JCl 



COT 

12A 

Xt'AI'IKtiXUM'Al}}! 



COT 

12S 

lei icMiAur 



EOT 

>2A 

NttyliKIV • l 

1 



CUT 

>IT 

1 

U FINE HI AtlA.nVB OINmv. MNinV. m Phsm VAAIAUOH bn tACtt 

cut 

12« 

1 HCAUfMAl IIHC 



CU 

>24 

re 1 U'«i,N4vi 



COT 

»0 

HFlHi' • t 



COT 

»t 

- rUM/FAEiUA 



cut 

m 

PC i A»l»NHP*t| 



CUT 

>» 

t riNG Ite A,¥.l CCOdCIAAIIS Of BACH «*0|KI 0N iMB UNB 



CUT 

»A 

FKitA 



CUT 

»s 

tt.ElAP»PlDIV? ♦ IfA ’> XHOaUM«tLN6/IH'K| 



CUT 

IIA 

cc$uNa.osnHptAPt 



CUT 

»T 

SlN»t‘P-BSlNM»<BTAt't 



CUT 

>)S 

KP»«*StMHl'»CCSll»F 



CUI 

>14 

yi'»A*StMHP»StMfP 



CUT 

2A0 

{p»A*cpsnip 



cut 

TAl 

A»AU*AP « Aiant' ♦ AliAlP 



cut 

IA2 

y«ASl*>iP * Ajjtyp * AH*tP 



CUI 

TAJ 

# AJ?»VP • Ali*lP 



cut 

lAA 

W IS THE I'OlNi lA »KC ElKBt 



CUT 

IAS 

ir lA " KtOAi) 



CUT 

>AA 

H JhUC»i2 



CUT 

lAT 

W 1 IS ESni FPA EPAf 



CUT 

lAA 

SC »C ill 



CUI 

>A4 

it CCNUNUe 



CUI 

ISO 

C ts llif: tMtlNI ftClOti AUAERt 



EOT 

>Sl 

ir lA - Aip^rm 



CUT 

l»2 

lei crN»o.iM‘o 



CUT 

)S> 

&r tc ic) 



CUI 

ISA 

ic? ifNi imit- 



CUT 

>» 

t riNC Ml SPtlBMCAl MUAA (OOAOlNATii OF POINT 



CUT 

ISA 

AS«TSE*CS0«1l**»it * ViAfl 



CUT 

1ST 

If lit iA,JS,SA 



CUT 

SSB 

J8 lttlA»(UDIV? 



CUT 

114 

C.C U it 



CUI 

lAQ 

H lltASEYSlVi 



CUT 

lAt 

Tt'BIA»lUTANtCT 



CUT 

IA2 

If lit AA^AiBA 



CUT 

lAI 

AA IPBIA.PI ♦ TIBIA 



CUI 

lAA 

Jt CCNIlNUe 



CUT 

IAS 

CAU ANGlBUiV.AlAOAl 



CUT 

>IA 

U FINE Tl't AtlATIVB lUNJIlV, OBNSITy. ON 0BN5ITV VAAIAIION 

AT POINI 

CUT 

lAT 

CAU CKASTV IA,1HtTAiAEAOAiNO«N,RrfOBN»OtNI 



CUT 

lAA 

ICI CCNIISUI 



CUT 

IA4 

E I'WT THE AlCiHT STNAOl/fUASlA IN AAAAIM/TNAAIM 



CUT 

JTQ 

AS«ICt.N/SCAEBI « C»!OC 



CUI 

in 

aMlAiSVAAE ♦ IS 



CUT 

IT2 

If IJMIA .E*. II JNlA-l 



CUT 

111 

If UKUA .C«. NAASYAI JAUN»AAISTA 



CUT 

ITS 

24 lAAAIAI'AEAl'AUAUAI 



CUT 

ITS 

TAAAIR) * ASTAtJNtNI 



CUT 

ITA 

1 ECNIINUB 



CUT 

ITT 

E I’AIM li t UNt 



CUT 

III 

NAtU lA.lAI IXPAI'IM, A*ltUU 



CUT 

IT4 

VA fCAAAV IM*12UU 



CUT 

lAO 

U NAIU IIO.lCCUIAAIMlUt K<I«UU 



CUT 

IKt 

t IQO rCAAAl I12U2I 



CUT 

>A2 

C AlStT ll-e UNB 



CUT 

a«i 

oc IS A«ia?i 



CUT 

lAA 

*PAl'IM«AEP|iAI2Jt 



CUT 

lei 

TAAPIAl • hhaivvi 



CUT 

lAA 

U CCNIINUI 



CUT 

JAT 

2 ECNT IKUt 



CUT 

in 

fei CGNIINUl 



CUT 

>«4 

AtTEAN 



CUT 

>40 

BKC 



CUT 

141 

SlSAACEtlNu Piece IAl.A2,Tlt8IAMI<eTA2*KEPnAl*XEPl)A*iNSPAf,R,N$llAFI 

e 

Pie 

1 

E 

E APieCE* 

r 



Ptl 

2 

t ti‘is sueAouuNt ccNfUies iiit AeiAtivi o«Nsm* dinsitv 

• DA. 

Pll 

1 

U CCNSI1V VAAUTIQN ON A fOAHIA Of SPHEAICAE AtCTANCEES* 

AEE 

PITH TUI 

Pit 

A 

e SANt ANCEEAA CCOAIUNAIES EOT POSIUCNtI) AT OlFWiNI OtPIllS. 

r 

Pie 

S 

E A SPflAlEAE AeCIANCU U Till AAtA ON THt SOATACB UP 

A 

Sfl'IAB 

Pit 

A 

E eOUNtHO BY INO rUtC EAll totes ANQ TAO FI ABO UtNOTIOOBS 

E 

l; SPHtAltAE HiCTANOltS AAt AENAYS COHPOTeO BOA AAOIAE 

‘ 


Pie 

1 

OISIANCeS At 

Pie 

( 


4a 


MO KIi hITK m OIHCRS (IF MVl SANCNICHEO IN lEmilN IMEJI fNO. 


Fli 


TNE SFHEMC4L KECIAfClEi ME CISPL4IE0 AS RECrANeiES 

WITH THE 

FIE 

10 

: EASI>«iEST OIAECTION RVNAIKC NOKIIONTAi|,V ACEOSS INE PACE 

IPFOTOI. SO 

PIE 

11 

INEAE IS A CEAIAIN AHCUAT Of DISTOKUON. 


PIE 

12 

NCTATIQN 


PIE 

U 

Al > AAOIAL DISIAKCt CP LOMEST IIECTANCEC IN PILCNETEAS 

PIE 

14 

At « AADIAL OISIAACI OF HtCHEST AECIANCLE IN KIIOAEIEAS 

PIE 

IS 

INEIAI • CCLAtItUCE OF AOATMEAN ECUMCAAV OF IHE AECTANClfS IN 

PIE 

14 

UEOAEES 


PIE 

IT 

TNEIA2 • COLAmUEE OF lOUTHEAN EQUKDAAV OF IME AECtANClES IN 

PIE 

II 

DE6ACES 


PIE 

10 

XLNCAI ■ iONCIFUOE CF IFE MES1EAK AIIUNUAAV OF IHE AECrAACLES IN 

PIE 

20 

DECAIES 


PIE 

21 

AINCA2 • LONCIIUDE CF FFE EASU'AN BDUNDAAY OF |»IE AECfAACI.ES IN 

PIE 

22 

OECAEES 


PIE 

23 

NSFACE • LENCIH OF IHE AOAFHEAN ANO SQUfHEAN AOUNOAAIES IN 

PIE 

24 

NOAIIONfAL SPACES ACAOSS HIE PACE IFHOfOI. 


PIE 

2$ 

NSUAF • NUPFEA OF FECI ANCLES OESIAltO 


PIE 

24 

SUOACOTINES CALLECI 


PIE 

2T 

1. CINSfV 


PIE 

28 

NCTES 


PIE 

29 

1. A> .Cf. Al. 


PIE 

30 

2. IPETA2 .Cl. fPElAt, 


PIE 

31 

1. ALNCA2 .Cf. XLPCAl. 


PIC 

32 

4. NSFACE .LE. 120. 


PIE 

31 

S. NSUAF .CT. 1. 


PIE 

34 

SAMPLE INPUf DAfA (COLUMN 1 OF INI'UI SIAAfS IN CCLUMN 3 HEAE) 

PIE 

35 

SUANME IFOkpATi aa) 


PIE 

34 

PIECE 


PIE 

3T 

ICOPENTIJIt J«lil4l (FOAMAfl llAA. 

A2) 

PIE 

38 

PANAIi 


PIE 

39 

Al A2 INETAI TIIEIA2 HLMOAI XLN0A2 

nspace NSUAFPIE 

40 

IFOAMAfl AFIO. 

S.2ISI 

PIE 

41 

.2121.0 N121.0 SC.C RO.O 180.0 2L0.0 

80 

3PIC 

42 

IMPLICIT AEAL*8IA>H,0-2I 


PIE 

43 

INTECEAP2 IMAPluili EEIAIIII 


PIE 

44 

DIMENSICN ALPNAI23),LCMFI(20S«2I 


PIE 

43 

OIMENSICN AMAPIUII 


PIE 

44 

DIMENSICN COPENTII4I 


PIE 

4f 

CCPPCN/ELAA/APAP 


PIE 

48 

CCPPCN/ELAC/F,P|.PI2.PlCIV2.$CALE.NIiHQrOiMAXS'ifM 


PIC 

49 

COPPCN/ELKDFALPHA.LCMFI 


PIC 

SO 

CCMPCNFBLNC/LPINfLMAX.LLOHtLUP 


PIE 

SI 

CCMPCN7RLAH/AADAVC.PCCAE,AL0HEA.AUPI>EA 


PIE 

S2 

CCPPCN/ELAI/CCMEN! 


PIE 

S3 

CCMPCNFELAJ/NOENtAEFOEN 


PIE 

S4 

CCMPCN/ELAf/AETA 


PIC 

55 

PAINT INPUT DATA 


PIE 

54 

NAME lAtSOOl 


PIE 

5T 

ICO FCAPAT l/FF«UX»'SUEACUriNE CALLEOi PIECE". ///) 


PIE 

58 

MAITE lA.ni) ICOPCNTI.)!. J«t.l4l 


PIC 

59 

)U FCAPAT (IOA.UAA.A2.///I 


PIE 

40 

FINC PICDLE OF STPACLACCLCA AANCE 


PIE 

61 

SVPPAX«PAXSYP 


PIE 

42 

$YPEL<>SYPMAN/2.00C • C.SOO 


PIE 

43 

STAECP*S.O0O/4.0OC 


PIE 

44 

IF INPHCTO .EQ. n STAECIW.OCO 


PIE 

4$ 

WAITE l«tl» 


PIE 

44 

1 FCAPAT tIHI) 


PIE 

67 

NSPt«NSPACE * 1 


PIE 

48 

SUAFCE-NSUAF 


PIC 

49 

SPACEMSPACE 


PIE 

TO 

OA«<A2 - A1>/(SUAFCF > 1.000) 


PIE 

71 

DLMCAMALHDA2 - XLMCAII FSPACE 


PIE 

72 

C1PE1A>CLM0A*STAECH 


PIE 

73 

NTHETA«ITHETA2 “ THETAII/OTHETA ♦ 0|.500000l00 


PIE 

74 






44 


ooors r»or&r»r>or»r>iT5nor?»r»ram 


Nim»Nth6l* * I 

* i»<»tneu 

C 0« * Clit*NCC DEIhlilN SU»U£t$ 

t PtHtT* » l*H!l,CC BgtWIth NtlCHBQRlNO PQINIS 

W KIRCA • LCNOlaCl dlHi:«fNl.e NEHJ»‘POHING POINfS 

*, 

hRMC (t)2) K|tK?tUH 
7 : FCPPAt UR. 'HI** iF»C.l( IS. VAP'.ICX. 

. nt.j.is.'sps//! 

MRttt U.)l UtCTAI. IM(U2.tnHCU 

I FCRPAI UX,MHeiAI»'.ri{*A.IS.'RECHfeS'.IOX.MHi:tAJ*».FIO.*.l*t 

. *DtCHCCS»*lOS. •DIHn*»*,riO.S,l»,‘U8CREKS*,//l 
ht>|tC It. SI XI.Hn»|,HPOU.DLPCA 

S rCHPAI UK.‘SlPC»l*'.flf.t.l»,*CICCRe(iS*,|OX,'ILPD»«*SFIO.A,lS, 

. »ctcptt$».ic«|tt!iPiA»'.no.A.is.*BraRi:E5*.//i 

I 

t NCHKAUIC OlSTANCtS U tAPIM RACIUS 

HI*Pl/HAOAVC 
R7I.II7/HAOAVC 
dR«CP/RA»AVC 

U eCNVEPI ASCl.ES fKCP OCCPEIS TO RAt'IANS 
i>ciAi»n.C'iAi»f 
TFE1AJ#!MEIA*»r 
im<nA«CTHEtAtF 
RLHCAl«iiL«OAI*r' 

*lHCA?»«l.Ht)A?»P 

OEPCA^OlHOAtr 

C 

WRITE It, II 

C 

U INinALIEE LINE VAI.I.E: 

C SWAP TS FCK PRIM-OU PICUIRO 
t TRAP IS FCR OICCREC PMHU 
CC A K*l.l2l 
rPAPIRItnCTAIll 1 
A *PAPIR1»AUI'PAUJI 

e 

R»RJ * C« 

I. 

CC 1C I SURF* I. USURP 
R«R - RM 

MitTAwUeiAl “ CTPETA 

oc 1 irpeiA*i,Niin 

!HETA»THETA ♦ CTHClt 
SlRCA«REHOAI - OIPD/ 

oc t R»l.NSPl 
XIRCA.XIROA ♦ CIRCA 

C eCRPUIC TR6 RiHAIlVE tCKSITV. OCNSITY. OR UENSITY VARIATICN At POI 
CAU CESSTV IH.1HEU.Xl.R»A.NQEN,HLrnCN,QCNl 
C PUT IRE RIGHT SVRPOl/RURMTR IN XRAI'IRI/TRAPIM 

xs>iccN/scAic) . c.toe 

JKUR»SVR»U » XS 
IF IJNUR .UC. I) JNIH»I 
IF UNUR .QE. RAXSYM JMIR-NAXSTR 
29 XRAPIKI»AU>t<AIJKUR| 

TRAPIKI « OETAIJNURI 
a CCNTIKUC 

C PRIM TPE LINE 

WRITE It.lAI IXRAPIX), X»l,|2l| 
lA FCRRAT ISX.12IAI) 

C WRITE I IQ.IOOmWAPIKl. R*I.UU 

C too FCRRAT lUlUl 

9 CCNTINUE 
WRITE It. II 

10 CCNTINUE 
RETURN 
ENC 

SUBRCUI INC PCUCS IHFN.AXGNAXtNRAC.NPOLEI 
»PCUE$* 

THIS SUPRllUtlNE CCHFUTES THE REUATIVE DENSITY. DENSITY. CR 
CEN5ITY VARIATION CN A SURFACE OF CHCSEN RADIUS U.SINC THE PCIAR 
CRTHCCRAPHIC PROJECTlCN. 


IT nSPUAYS THE NCR til PtLAR REOICN 
PCUAR REGION tr* NPfiie»2. 


NCTATIQN 


IF NPCLCol and the SCUTII 


RRR • RADIAU OISiANCl OF RAP SURFACE IN RIUONETERS 
ANCRAX » ANCUUAR Cl STANCE FROR PCUF TC HAP PERIHEIER IN OCGRtES 
NRAC « RADIUS OF THE RAF IN HORUQNTAE SPACES ACROSS THE FAGF 
iPHCiqi 

NPQIE • I FCR NDRTI FOlAR RAPi »2 FOR SCUTH PCUAr pap 


PIE 

TS 

PIE 

7t 

PIE 

IT 

PIE 

T8 

PIE 

TV 

PIE 

«0 

PIC 

as 

PIE 

R2 

PIE 

RT 

PIE 

RA 

PIE 

as 

PIE 

at 

PIC 

at 

PIE 

aa 

PIC 

av 

PIE 

VC 

PIE 

VI 

PIE 

V2 

PIE 

VA 

PIE 

VA 

PIE 

95 

PIE 

Vt 

PIE 

VT 

PIE 

VB 

PIE 

VV 

PIE 

ICC 

PIC 

loi 

PIC 

102 

PIE 

10) 

PIE 

IQA 

PIE 

105 

PIE 

tot 

PIC 

lOT 

PIE 

loe 

PIC 

lOV 

PIE 

lie 

PIC 

Ul 

PIC 

112 

PIE 

lU 

PIE 

IIA 

PIE 

115 

PIE 

lit 

TPIE 

lU 

PIE 

u« 

PIE 

uv 

PIE 

120 

Pl.E 

I.TI 

PIE 

12T 

Pl« 

12) 

Fi.e 

12A 

PIE 

125 

PIC 

126 

PIC 

UT 

Pie 

laa 

PIC 

uv 

PIE 

UiO 

PIC 

m 

PIE 

uu 

PIC 

III) 

PIE 

t'lA 

PIE 

II 5 

PIE 

Ut 

POL 

1 

PCU 

2 

POL 

) 

PCU 

A 

POL 

5 

POU 

t 

POU 

7 

POU 

a 

POU 

V 

POU 

to 

POU 

u 

POU 

12 

POU 

U 


45 



c 


C SURRCUriNCS C/iUtCt 


POL 

14 

1. 1. «hClE 


POL 

IS 

C 2. CIRC 


POL 

14 

C 3. CCM5TV 

C 


PCI 

17 

c 

c Ncris 


POL 

II 

Q 

C 1. ANCMX .11. 9C CECRCES. 


POL 

19 

C 2. .LE. «C. 

C 


POL 

20 

C 

C S2HELE INPUT D*T* (COtUHN 1 OP INPUT STARTS IN 

g 

CCLUPN 1 HERE) POL 

21 

C SUBNHE IFOANATl 

A6| 

POL 

22 

C POLES 


POL 

23 

C (CCPENTIJI. J'iflAI (FOAPATI 

I3A6.A2) 

POL 

24 

C SURFACE CENSITV CISTNleLTlOF 


POL 

2S 

C NRAO ARCNAR RNP IFORPATl 

IS.3F10.M 

NOL 

26 

C 3A 9C.0 tJTI.C 

C 

r 


P05 

27 

V 

c 

IPPLICIT REAl*e(A«H«0>n 


POL 

21 

INTECERP2 TPAP|t2ll« EETAUlt 


POL 

29 

OIPENSICN ALPHA(23) |LCWFl|20St2» 


POL 

30 

DIPENSICN MPAPU2U 


POL 

31 

oirensicn copentiiai 


POL 

32 

CCPPCN/ELKA/XPAP 


POL 

33 

CCPPCN/FLRC/FiPI.P|2.PICiy2iSCAL£.NPHOTO.NAKSVP 


POL 

34 

C,CPPCN/eLKD/ALPHA,LthM 


POL 

35 

CCPPCN/ElRC/lPIN,lP2)tfUOti,LUP 


POL 

36 

CCPPCN/EIKH/RADAVC»FCCRE>RLOWERiRU.PPER 


POL 

37 

CCPPCN/ELRI/CCRENT 


POL 

3B 

CCPPCN/eLKi/NOENtREFCEN 


POL 

39 

CCPPCN/eiRT/PETA 


POL 

40 

e PRINT INPUT DATA 


POL 

41 

WRITE Utl6! 


POL 

42 

WRITE (6.100) 


PCL 

43 

ICO FCRPAT (///.lOR.'SUERCUIINE CALLED.* POLES'.///) 


POL 

44 

WRITE 16. Ill) ICOPENTlJIt J>1.IA) 


POL 

4S 

111 FCRPAT I10X.11A6.A2.///I 


POL 

46 

WRITE 16.301) RRM 


POL 

47 

301 FCRPAT (lOX.'RKP •' .1*. FIO. 3 , IX. 'KR* . lOX »• (RAO 1 AL DISTANCE 

OF SURFPOL 

41 

.ACE)'.//) 


POL 

49 

WRITE (6.302) NRAC 


POL 

SO 

302 FCRPAT (lOX.'NRAO • '. IX .1 4 . lOX. • (RADIUS OF PAP IN H0RI2CNTAL SPACEPOL 

51 

.SI'.//) 


POL 

52 

WRITE 16.303) ANCPAA 


POL 

53 

301 FCRPAT (lOX.'ANCPAX IX. FIO.S. IX. 'DECREES'. lOX 

:.' (PAxiPUP 

COLATITPOL 

54 

,UCE)'.//I 


POL 

55 

R>RKP/RAOAVC 


POL 

56 

ANCP'ANGPAXRF 


POL 

57 

SANC'CSINIANCPI 


POL 

51 

C FINC PICDLE OF SYPECL/CCLCR RANGE 


POL 

59 

SVPPAX.PAXSyP 


POL 

60 

SYPeL*SVMMAX/2.00C * C.!DO 


POL 

61 

C CCPPUTE NUPBER OF HCRI2CNIAL LINES IN PAP 


POL 

62 

NRC5£C»NRA0 


POL 

63 

RACSEC*NRAD 


POL 

64 

NR¥»A»NRDSEC/5 


POL 

65 

IF (NPHCTO .EB. 1) NRT>NROSEC 


POL 

66 

NY«2RNRY ♦ 1 


POL 

67 

NRYloNRY « 1 


POL 

61 

••NRYI 


POL 

69 

CC • ..«c LEFT AND RIGHT LIPIT! OF PAP BY CALLING CIRC 


POL 

70 

CALL CIRC(NRCSEC) 


POL 

71 

DC 10 NUPPAP'I.I 


POL 

72 

Q 

C INITIALI2E LINE VALLE! 


POL 

73 

C XPAP IS FCR PRINT-OLT PICTURE 


PCL 

74 

C TPAP IS FCR OICOPEO PFOlO 


POL 

75 

DC A K*I.121 


POL 

76 

TPAP(X) • BETA(II) 


POL 

77 

4 XPAP|X).ALPHA(23) 


POL 

78 

Q 

WRITE (6.161 


POL 

79 

16 FCRPAT (IHI) 


POL 

10 

IF INPOLE .EB. 1) WRITE (6.30) 


POL 

11 


46 


r>o ft o r» ft ft r* ft r»ft ftftftoftooftftftftftft or* 


ir (K)>gu .ce. it hi>iie UiHI 

3 C IIOKf'NOMM ( 0 U»,/////i 

H rquMt UPXi'SOVTH roiE'./////i 

; EINC lEE «El«fiVE OtNSnVt CENSIIV, OR DENSI TV VARIRIIOM CN 
C HCRUCKtRL LINE 
OC J IR-l.NY 
LlpRV -If * > 

K>>IR>LOb> 4 mi*,n 
KRARpLOkHIUPiO 
DC 3 KpHRINiRHRX 

; MNC fh( M,yu CaCRC|RRU$ OR E«CI* POtNf ON THE MNE 
ERpR 

R«FK > euocc 
STRtCN*S. 000 / 4 . DOC 
IF INFHCIQ .CO. II ilREOlPl.OOO 
V«im p RY|l*StRICF 

: FINC Iht SPHERICAL FOUR (OOROlNAtES QF FO|Nr 

CALL ANCLCm.YfldFD/) 

IF INFQLE .EC. 21 XLHCA«P |2 - XLNOA 
X$CYSC>CS 0 RI(X «*2 ♦ Y»»JI 
SlNtHE>X$CY$G*$ANC/FAC$tC 
IF ISINIHE > t.ODCI 2 C 2 « 20 li 201 
HRIIE 14.2031 NRAC.FRCSIC.NY.NRY.NRYI.LF.L/.K 
FCRFAT I 3 X. III 3 I 

HRIIE U. 204 I RAOSet.RYl. FL 2 .fr. A. Y.X 50 YS 0 . SANG. SINIHE 
FCRRAT ISK.IIFIO.!! 

CC 1C to 
CCNIINUE 
CCSIHC»CSORTIl.ODC - SIMHE»»2I 
CALL ANCLEICOSIHE.SlNIHE.tHETAI 
IF INFULE *£«. 21 TFEIA.PI - IHEIA 
C CCHFUtC THE NELAriVE lERSITY. DENSITY. OR DENSITY VARIATICN 
CALL CERSTY IR. THETA. ALFDA.NDEN.REFDEN. DENI 
C PUT THE RIGHT SYPFat/FUPOIR IN APAPIR I /THAPIR | 
X5«ICEN/SCALEI ♦ C.fDC 
JNUPpSYPOL • *s 
IF IvINlIP .LE. II JNLP*1 
IF IJNUP ,GE. NAA5YM JRUK-PARSYP 
24 XPAPIKI>ALHHAUNUPI 
TPAPIRI • RETA(JNLM| 

3 CCNTINUE 


EACH 


201 

2C1 

204 

202 


: PRINT THE LINE 

RRITE It.l4» IRPAPIR), R»l,|2l| 
14 FCRPAT tSR.IOIAtl 
; WRITE IIO.IOCHTNAPIRI, R>I.I21| 
: too FCRPAT II2II7I 


13 


RESET THE LINE 
DC 13 R-RHIN.RHAX 
RPAPIRUALPHAI23) 
TPAPIRI - OETAIIII 
CCNT INUE 


2 CCNTINUE 

"iC CCNTINUE 
RETURN 
ENC 

SUnPCUTtNE SECTNSINSEC.NRAO) 


•SECTNS* 


THIS SUBROUTINE CCPfUTES THE RELATIVE DENSITY. DENSITY. OR 
CENSITY VARlATtON CN A NUPBER OF PLANES WHICH SLICE THRCUOH THE EARTH 
PARALLEL TO TILS EGUATCR. THE INTERSECTION QF THE EARTH WITH A PLANE 
IS CALLEC A SECTION. 

ONE SECTION ALWAYS CONTAINS THE EQUATOR. THE SECTICNS ARE EOUALLY 
SPACED BETWEEN THE NORTH POLE AND THE SOUTH POLE. FOR ERAPPIE. IF 
THREE SECTFCNS ARE CHOStN. CNE CONTAINS THE EQUATORIAL PLANE AND THE 
OTHER TWC SECTIONS ARE IPACCO ONE-HALF EARTH RADIUS NORTH AND SOUTH 
CF THE EGLAIOR lASSUPINt NO CRUST IS STRtPPEC OFFI. 


notation 


NSEC 

NRAC 


NUPBER OF SECT IONS CHOSLN 
RADIUS OF iSLATORIAL SECT! 
PACE IPHCTtI 


SUBRCUTInES CALLECt 

I. ancle 


POL 

«2 

PCL 

13 

PCL 

14 

POL 

«3 

PCL 

•4 

POL 

• T 

PCL 

■ 1 

POL 

•4 

POL 

40 

POL 

41 

POL 

42 

POL 

43 

POL 

44 

POL 

45 

POL 

44 

POL 

4T 

PCL 

44 

POL 

44 

POL 

too 

POL 

toi 

POL 

102 

POL 

103 

POL 

104 

PCL 

lOS 

POL 

104 

POL 

lOT 

PCL 

tOB 

POL 

104 

PCL 

110 

POL 

III 

POL 

.112 

PCL 

113 

TPCL 

114 

PI!L 

US 

POL 

l|4 

PCL 

in 

PCL 

Its 

POL 

114 

POL 

120 

PCL 

m 

POL 

122 

PCL 

m 

POL 

124 

POL 

123 

POL 

124 

PCL 

12T 

POL 

I2B 

POL 

124 

POL 

130 

POL 

III 

POL 

112 

POL 

113 

POL 

114 

POL 

113 

POL 

114 

POL 

IIT 

EEC 

1 

SEC 

2 

SEC 

3 

SEC 

4 

SEC 

3 

SEC 

4 

SEC 

7 

SEC 

8 

SEC 

4 

SEC 

10 

SEC 

II 

SEC 

12 

SEC 

13 

sec 

14 

SEC 

15 

SEC 

14 

SEC 

IT 


47 


o rsr»r^or>oc^r>o<ir»r»r" nr^ or*’ r>oc*<^ 


*. CIRC sec ve 

J. ccNstv sec IS 


NCTES 


SIC 

20 

t. Nsec 

IS ccc. 

SEC 

21 

if RR*C 

.ie» »c« 

sec 

22 


StKIU INPUT C*TA ICCUHN ! 

1 01 INPUT STARTS IN COUNN 3 HEREI 

SEC 

23 

SUHNRI 

IIORHATt A 

SEC 

2A 

5ECTN5 


SIC 

2S 

ICORiNTIJt, 

IFORHAIt 13AA»A2I 

sec 

26 

live SECTIONS ChCSEN 


SEC 

27 

NSEC NRAC 

IIORPATl AISI 

SEC 

2« 

S 6C 


sec 

2« 


IPPLICIT REAIA«(*-M,0>n 

SEC 

30 

INTtCtR*2 TPAPrUll, EEIAUH 

SEC 

31 

OIPENSICN AlPhA|23liLCWHI205.2l 

SEC 

32 

CINENStCN XPAIM12I) 

SEC 

33 

DIPENStCN CCPENTI lAI 

SEC 

3A 

CCPPCN/eiKA/APAP 

SliC 

35 

CCPPCN/Pt.>!C/F.PI>Pi:.PICIV2.SCALe.NPH0T0,MAJlSyH 

SEC 

36 

CCPPCN/eCKD/ ALPHA, UNH 

SEC 

37 

CCPFCN/DLKC/LPlK,LP/K,LLOH,LUP 

SEC 

38 

CCPPCN/eLK»l/RAOAyC,FCCRE, RICHER, RUPPEH 

SEC 

30 

CCPPCN/eLKI/CCPENI 

SEC 

AC 

ccppcn/dlkj/ncen.retoen 

SEC 

A1 

CCPPCN/eLKT/BETA 

SEC 

A2 

: PRINT INPUT CAIA 

SEC 

A3 

NRITE U,16I 

SEC 

AA 

WRIIE 11,3001 

SEC 

AS 

ICO FCRPAT l///,VOX,<SUIRCUUNE CAUEOi SECTNS',///! 

SEC 

A6 

WRITE 16, Mil ICOPEMIJI, J*l,lAI 

SEC 

A7 

311 FCRPAT llOX,nA6,A2.///l 

SEC 

A8 

WRITE 16,3011 NSEC 

SEC 

AS 

3Cl FCRPAT I10X,'NSEG '•MJt,lA,IOR,'|NUIiDER CF SECTICNSI* ,//) 

SEC 

5C 

WRITE 16,1021 NHAC 

SEC 

SI 

302 FCRPAT llOit,'NRAO • », l« ,| A, lOR, • IRAD lUS Cl EOUATCR IN HCRUQNTAL 

SSEC 

52 

.PACESI',//) 

SEC 

S3 

: riNc PicoLE OF svpbcl/ccur range 

SEC 

SA 

SyPPA*«PAi(SYP 

SEC 

85 

SVPei«SyPNAX72.0DC « C,100 

SEC 

56 

NHEPIolNSeC - l)/2 ♦ 1 

SEC 

57 

FRAC'NRAn 

SEC 

58 

MEPIoMIEWI 

SEC 

51 

i DC EACH SECTION 

SEC 

60 

OC 1C JJ«l,NSEC 

SEC 

61 

FJJ-NHEPI - JJ 

SEC 

62 

IF IJJ (,E0, NHEHtI CO TC 21 

SEC 

63 

SCUARI«IF«A0*»2I»I1.0C0 - tFJJ»»2)/lHENI»»2| 1 

SEC 

6A 

KAOSECoCSaRTISCUAREI 

SEC 

65 

2R>iFRA0A|PJJ/HEPI 1 

SEC 

66 

NRCSEC-RAOSeC 

SEC 

67 

RADSeC'NROSEC 

SEC 

68 

CC TC 22 

SEC 

6S 

21 NRCSEC>NRAO 

SEC 

70 

THETA-PI0IV2 

SEC 

71 

HACSEC'FRAD 

SEC 

72 

22 CCNTINUE 

SEC 

73 

: ccwpuie nupoer of hcriecnial lines in hap 

SEC 

7A 

NRy«A*NR0SEC/5 

SEC 

75 

IF INPHCTQ .E0> 11 NRV>NRDSEC 

SEC 

76 

NV«2ANRy ♦ 1 

SEC 

77 

NRVI»NRy ♦ 1 

SEC 

7* 

RV1»NRV1 

SEC 

79 

FIND LEFT AND RIGHT UH|T< OF HAP By CALLING CIRC 

SEC 

80 

CALL CIRCINRDSEC) 

SEC 

81 

: INITIALI2E LINE VALUE! 

SEC 

82 

RHAP IS FCR PRINT-OUT PICTURE 

SEC 

13 

: THAP IS FCR OICCHEO PFOTO 

SEC 

8A 

OC A K«l,12l 

SEC 

85 

TPAPIKI • BETAUII 

SEC 

86 

A XPAP|«|*ALPHA|23I 

SEC 

87 

: SPACE CCWN THE PACE 

SEC 

88 




♦ 


48 



UHltE lEiltt 



SEC 

•5 

u 

PCHPAI UHU 



sec 

90 


NSPACe*l2« «■ NftV 



SEC 

91 


DC 17 KK«l.NSPACE 



SEC 

92 

n 

Num ltil03 ALPhAUJ) 



SEC 

93 

le 

rropAt iioXfiAU 



SEC 

99 

V 

2.PjJ»«LPP£R/HEPt 



SEC 

95 

c 

FINE TPE KELAIIVE DfN!l!V. DENSlIVi ON DENSITY VANIAMON 

CN EACH 

SEC 

96 

L 

HCRMENTAl, LINE 



SEC 

97 


DC 2 LP-l.NV 



SEC 

9« 

C 

FINE TFE YiY«2 CCCRCINA1E! QF EACH POINT 

ON the line 


sec 

99 


U«NY -LP • 1 



SEC 

too 


FU»Lt 



sec 

tot 


STRECP»5.0D0/A.CCC 



SEC 

102 


IF IKPHCrU .ED. It !TNECH*l.0C0 



SEC 

t03 


Y»iril - RYll«l5l»CtHARlPPE»l/RACSEC 



SEC 

109 


KPIMLCNHIIlPill 



SEC 

tos 


KPA«>Lbt>MIILP«2l 



sec 

t06 


DC 3 KaRHINiKPAX 



SEC 

107 


F*»N 



sec 

to« 


*a|FN - 6l.qC0l*|RUfPER2HA05ECI 



sec 

t09 

c 

FIND TFE SPHERICAL POLAR CCOROINAIES OF 

POINT 


sec 

uo 


R»DSCRnx»«2 ♦ Y*«2 ♦ 2*»21 



SEC 

III 

c 

IS THE POINT IK THE CCRE7 



SEC 

112 


IF IN ■> RCDRCI 2Ef27(27 



SEC 

113 

2e 

aNUP«22 



SEC 

U9 

c 

» IS liSEO FCR CCRE 



SEC 

US 


CC TC 29 



sec 

116 

27 

CCNTISUE 



SEC 

117 

C 

IS THE POINT RELOH RLCHIR1 



SEC 

III 


IF IN > NLONEHI lCl.tC2.t02 



SEC 

119 

iCl 

DENaC.DCO 



SEC 

120 


CC 1C IC3 



SEC 

121 

IC2 

: CCNTINDE 



SEC 

122 


XSCYSC«CSgRTI***2 ♦ YA*2I 



SEC 

123 


IF HI 2A. 25,29 



SEC 

129 

25 

TPETA«p|0|V2 



SEC 

12$ 


CC TC 26 



SEC 

126 

2^ 

U«*SCYSC/2 



SEC 

127 


THETA-OATANILi 



SEC 

12R 


IF m AA.2At26 



SPC 

129 


THETAaPI ♦ IFEiA 



sic 

130 

26 

CCN1INUE 



SEC 

131 


CALL ANCLnX.Y.XLPOAl 



SEC 

132 

C 

CCPPUIE TPE RELATIVE lENSITY. DENSITY, OR DENSITY VARIATICN AT POINTSEC 

133 


CALL CENSTY (R , THETA, XL PDA .NOEN.REFDEN 

.DENI 


SEC 

139 

IC3 

CCN1 INUE 



SEC 

13S 

C 

PUT TH RIGHT SYNPCL/PUPCER IN XHAPI K) /TMAPIK 1 


SEC 

136 


XS*(CEN/SCALEI * C.50C 



SEC 

137 


JNUP*SYPHL • XS 



SEC 

131 


IF IJNUP .LE. 11 JNLP«l 



SEC 

139 


IF IJNUP .CE. PAXSVPI JNUPaPAXSYP 



SEC 

190 

29 

XPAFIXIaALPHAIJNUPI 



SEC 

t91 


TPAPINI • BETAIJNLPI 



SEC 

192 

3 

CONTINUE 



SEC 

193 

C 

PRINT THE LINE 



SEC 

199 


WRITE I6,IA| IXPAPIPI, P-l,l2M 



SEC 

195 

u 

FCRPAT l5X,t21Al) 



SEC 

196 


HRITEUC.tOOUTPAPiP), P«l,l2t) 



SEC 

197 

too FCR«>M U 2 U 2 ) 

r 



SEC 

196 

c 

RESET THE LINE 



SEC 

199 


CC 15 RaKPIN.KPAX 



SEC 

150 


XPAPIRI.ALPHAI231 



SEC 

151 


TPAPINI « RETAtlll 



SEC 

152 

15 

r 

CONTINUE 



SEC 

153 

U 

2 

CONTINUE 



SEC 

159 


XLAI»PHI«lt80.0CO/Pil 



SEC 

155 

C 

PRINT LATITUDE DF EACH SECTSON 



SEC 

156 


WRITE (e,30| XLAT 



SEC 

157 

3C 

FCRPAT l///,90K,'LAmUtE>',F9.4,IX.'DEGREESM 


SEC 

151 

tc 

CONTINUE 



SEC 

159 


RETURN 



SEC 

160 


ENO 



SEC 

161 

c 

SUGRCUTINE SLICE IPHI 1, PLPDA 1 ,PKI 2 , XLPDA2, TOLERO.NRAD 1 


SLI 

1 

c 

c 

•SLICE* 



SLI 

2 

c 

THIS SUBROUTINE CCPFUTES THE RELATIVE 

DENSITY/DENSITYAOCNSITY 

SLI 

3 

C V*RIATKK ON A PLANE hH ICP PASSES TNNOUCH 

r 

THE CENTER OF THE 

EARTH. 

SLI 

A 

c 

/* 

THE INTlRSECriCN CF THg PLANE WITH THE 

EARTH IS CALUEO 

A SLICE. 

SLI 

9 

c 

THE SLICE IS ORIENTED SC AS TO CCNTAtN 

THE TWO PCINTS 


SLI 

6 

C (PHIt.XLPCAU ANC (PHI2,XLNCA2)ii IHE LliNE 

WHICH FORMS THE 


SLI 

7 


oonooor»r>or»or>norvr)oor>or>ooor»/^or»rw*>r^rior»oofvor>or>ir»oinor>f> 


INrCRSCCIlOK OF fNE SLICE WIIH THE tCUAIOR RLNS HORItONTALLt ACROSS 
THE FACE IFHOTOI. NE ARE lOCRINC DOhR ON THE SLICE FROF THE RORTHERH 
HEHISFHERE* 


ROTATION 

FHII « LAlirUGE OE FIRST ROIM IN OSGREES 

KLRCAI • LONCIIOOE CF FIRST ROINf IN DECREES 

RHI2 « LATITUCE OF SECOND ROINf IK DECREES 

XLRCAE • LONCITLOE CF SECOND ROIRT IN DECRIES 

TCLERC > rOLERARCE lIRIf IN OECREESf EELCH NHICH A SHALL 

difference |R latitude CR LONCITUOE OF THE TNO ROINTS 
IS SET ECUAL TC EERO TO AVOID STNCULARITIES 
NRAC • RADIUS OF 5LICE IN HORUONTAL 5RACES ACROSS THE RACE 
IRHCTOi 


SUERCUTINES CALLECI 

I. ANCLE 
i. CIRC 
3. CENSTV 
A. VCC 


NCTES 

1. RLRCA2 .CT. RLHCAI. 

2. RRAC .LE. 60. 

3. SUBROUTINE CUTRLTS CAROS FOR CALCCHP RROCRAH. 


SANPLE INPUT DATA ICOLUPN 1 OF TNPUT STARTS IN CCLUHN 3 HERE) 

SUBNHE IFORHATl A6I 

SLICE 

ICORENTIJ). U>I.1AI IFORHATl 13A6.A2) 

SLICE RUNS CONN THE HAUAIIAR CHAIN FROM MIOHAY 10 BIG tSLANC 

PHIL XLHDAl PHI2 XLHCA2 TOLERO NRAC 

IFORHATl 5FI0.5,I5» 

26.11 163. 3S LS.cA 2C6.S9 C.COCl 6C 


IHPLICIT REALASIA-H,0'2I 
INTECER*2 THAPI12I). EETAIUI 
OIHENSICN ALPHAI23).LCUFU2CS,2) 

OIHENSICN XHAPII2I) 

OIHENSICN CCHENTIUI 
CCHHCN/BLKAFXHAP 

CCHRCN/eLKC/F.PliR[<.PIt:;V2.SCALEiNPHOTO,HAXSYH 
CCHRCN/BLKO/ALPHAiLCHHI 
CCHRCNFBLKC/LHINtLHAXtLlOH.LUP 
CCHRCN/ELKH/RAOAVCiRCCRT.RLChER«RUPPER 
CCHHCN/6LRI/C0HENT 
CCHRCN/6LKJ/NOEN.RETOEN 
COHRCN/BLRT/BETA 

C PRINT INPLT DATA 
NRITE 16.16) 

)IRITE 16.300) 

3C0 FCRRAT I///. lOX . • SUFRCU1INE CALLED! SLICE'.///) 

WRITE 16.311) ICORENTTJI. J«l.U) 

311 FCRRAT I10X.I3A6.A2.///) 

WRITE 16,67) PHII.XIHCAI 

AT FCRRAT I3X, 'FIRST GREAT CIRCLE PCINT I ' , SX, ' LAT I TUOE>' . lX.F LO.S. 
. IX, 'DECREES', lOX, 'LONG ITUUE»',1X,F10. 5, IX, 'DECREES',//) 

WRITE I6.A6I PHt2,XlHCA2 

A6 FCRRAT ISX.'SECCNC CHEAT CIRCLE POINT I ' , SX, 'LATI TUOE>' , IX ,F10.S , 
. IX,'CECREES',10X,'10NCITU0E>‘, IXiFIO.S, IX, 'DECREES' ,// ) 

C CCNVERT ANCLES FRCK OtCREES TO RADIANS 
PMI»PHIl»F 
PHI2"PH12*F 
XLHEAl»»lHOAJ*f 
XLPCA2<XIHOA2*F 
TCLER«TCL€RO»F 

C FINC XXI, ONEGA, COSXI. SINXt, COSCHC. SINCHC 
TCLPHI*C50BTIPHIl»«i • FHI2»»2I 
C IS THE CHEAT CIRCLE TFE ECUATGR7 
IF IICLPHI - TOLER) 7S,60,80 

/9 XXl'C.OCO 
CCSXI'UODO 
$INXt«O.ODO 
OHECA'O.ODO 
CCSCRC-l.OOO 


SLI 

1 

SLI 

S 

SLI 

10 

SLI 

it 

SLI 

12 

SLI 

13 

SLI 

lA 

SLI 

IS 

SLI 

u 

SLI 

17 

SLI 

11 

SLI 

IS 

SLI 

20 

SLI 

21 

SLI 

22 

SLI 

23 

SLI 

2A 

SLI 

25 

SLI 

26 

SLI 

27 

SLI 

26 

SLI 

26 

SLI 

30 

SLI 

31 

SLI 

32 

SLI 

33 

SLI 

lA 

SLI 

35 

SLI 

36 

SLI 

37 


SLI 

36 

SLI 

39 

SLI 

AO 

SLI 

A1 

SLI 

A2 

SLI 

A3 

SLI 

AA 

SLI 

AS 

SLI 

A6 

SLI 

AT 

SLI 

AB 

SLI 

A9 

SLI 

SC 

SLI 

51 

SLI 

52 

SLI 

S3 

SLI 

56 

SLI 

55 

SLI 

56 

SLI 

57 

SLI 

58 

SLI 

59 

SLI 

60 

SLI 

61 

SLI 

62 

SLI 

63 

SLI 

6A 

SLI 

65 

SLI 

66 

SLI 

67 

SLI 

66 

SLI 

69 

SLI 

TO 

SLI 

71 

SLI 

72 

SLI 

73 

SLI 

TA 

SLI 

75 

SLI 

76 

SLI 

77 


50 


SINO'CvC.OCO SU 

CC 1C 91 SLI 

HC rcu**'»cA#sn«u<c9« * *u 

i IS iMt cm*? ciHcu * HcmoMN) slj 

I* ncu*?' - ?Cl£*» 9l.fl.9J 511 

82 mimiciva su 

CCSM'O.OOO 51 

siN*i*t.oi)o stJ 

DK£C*»KlMU*l 5U 

cc$c*'e«ccosa*Ec*i sl» 

5lhC*0*CSIMtC»'EC*) *«•' 

CC 1C «l SC| 

19 CCNtlKUl , S1.I 

flNC*'C«UD5IMXU'CAII lilOTANinMIin * IOSIMiaf‘0*2l l*IOr*NIPHI I ) 1 1 SU 
./(ICCCSUamAin^lDUMPHiatt - !0COSr*l,l'DA2H*<01*NIPHUH} su 

OPfCA'CArANITANCPCI SLI 

IlNlXt»tDr*M(PNI7ll/(C5tNtXl.HC*2 • OKECAII 5ll 

If lUNXXII 92.52.5} SU 

52 OPtC»»0*EO* ♦ MI SU 

1*N)i»l»-r*NXXI SU 

53 CCNItXUC SCI 

X«l>C«r*NIUKXX|) SU 

IF (CPEC*I 55. 5*. 5* SU 

55 0*eC**0MEC* * P|2 SU 

54 CCMlFiUE SU 

XXIC*XXI/F su 

CPECAC^CMEOA/f SU 

: l»fUM XXI *N(J CMEO* IF CECREES SU 

: XXI « xncle of o«e*i u«cu mic ecuator su 

: OPEC* « tCNciacE cf /scefdinc acce of creaf circle. PEAHREn su 

: EASIRARC FRCP CREERWIOl SU 

WRIIE It.jCl) XXIC su 

JCI FCRPAI IIOX.'XKI ...ix.f 10. 5. IX. ‘DEGREES'. lOX.MINCLINAIlON TO EOUSU 

.AlCR)'.//) SU 

HRIIE IE. 1011 CPECAC SU 

103 FCRPAT IIOX.'QPECA •• ,1 X.F 10.5. IX. 'DECREES lOX. ' (NODE FEASUREO EASLI 

.SI FRCP CREERwICHI' <f/l SLI 

: PRIM ARAC SU 

WRITE IE. 1021 NRAC SLI 

3C2 FCRPAT IIOX.'NRAO « • . IX , | A , lOX. • IRAOlUS CF SLICE IP HCRIZONTAL 5PASLI 

.CESI'.//I SLI 

CCSXIxDCOSIXXIl SU 

$INXI>OSINIXXn SLI 

CCSCPC>CCOSICPECA) SLI 

SIN'CPC»CSINICPECA) SLI 

Bl CCNIlPUE SLI 


: CCPPL1E CALCCPP PLOIIER DATA FOR PIDITINC TOPOORAPHT OF SLICE SU 

C ARC GREAT CIRCLE CM VAN OER CRIN1EN PAP SLI 

C RACPAP > RACIUS CF VAN tER CRINTEN NAP IN INCHES SU 

RACPAP*I0.|RT500 SLI 

S CCNVER • INCH TC XN RATIO USED IN lOPCCRAPHY PLOT SU 

CCNVER'O.lOO SLI 

RAD5EC-NRAO SLI 

RCINCH-RAOSECPIC.IDCI SU 

RINCH'RCINCH » I.COE SLI 

NCEC«1 SLI 

WRITE IE.20BI SLI 

2Ca FCRPAT I//////. lOX. 'CALCOPP PLOT DATA',//) SLI 

WRITE IE.20AI RAOPAF SU 

2CA FCRPAT ||2X.'RADILS OF VAN OER CRINTEN NAP >'.F7.3. IX. • INCHES' ) SU 

WRITE IE. 2051 CCNVEF SLI 

205 FCRPAT I //• 1 2X. ' INCFE S PER KP Of TOPO HEIGHT ••.ne.51 SLI 

WRITE IE. 161 SLI 

WRIIE (E.20EI SLI 

206 FCRPAT I2X,'K'.7X.'FSI'.TX. 'LAIITUOE'. AX. 'LONGITUDE'. 2X. .TOPO HEICSLI 

.HI'.3X,.XGRIO'.7X,'TGR|E'.IOX.'XH'.IOX.'VH' I SLI 

WRITE 16.2071 SU 

2C7 FCRPAT |BX. ' lOECREESI ' . 3X. ' IDECREES I ' . 3X. • lOECREES I ' .5X . ' I KN I • . SU 

. EX.'MNCHESI'.AX.' IINCFESI'.5*.MINCHESI'.AX,MINCHESI'./I SLI 

OC 201 X>1.360 SU 

PSIC-R - I SU 

PST«PSIC*F SLI 

ccspsi»cca$4psn su 

SINPSI>CSINlPSII su 

X«CCSPSI*COS0PG - CCSXI*SINOP0»SINPSI su 

y-CCSPSIESINCPC * CCSXI‘'COSOPC*SlNPSI SLI 

E"SINXI*SINPSI su 

xsoysc«csuRTix»»2 ♦ y»«ii su 

IF IXSOySQI 66.67.66 SLI 

67 PHIC«PICIV2/F SLI 

THETA'O.ODO SLI 

IF III BA.89.85 SLI 

BA PHIC»-PI0IV2/F SLI 

THETA-PI SLI 


51 


ii 

i 


i 

i 


e$ cct4ll^uf 
CC TC i>« 

6t U«{/»SQV$0 

PKt»C»r*N«U» 

tKUi>PI01V? • PHI 
PHC*PHI/F 
6B CCMIWt 

CPU PNCce u.ytXLHtAi 
KIPC*C»»U»B*/F 
IF (ALHCADi ntn,n 
7J KlPtPC«»LHC*C ♦ Ho.qco 
7A CCNIlNUt 

CPU lOPOOA IIHCIPtHPCPiH) 

HF»I*CCAVE« 

CPU VDC tPHICiXUDIOtHIDHPPiNOeCUCiYC) 

XH-IPINCH ♦ FHI»CCSPSI 
VH.IBINCH * FHI'SINfSI 

2C2 rCKPPr llKtlPi»l2PinC.2MP.FlO.)i2XiPtC.1l 
HPHE le»202) KtP$ll«PHIO»PlPC*DtN.KC»YC(PH«YH 
UPltC 11)2011 I()XC)YC)PF)YM 
2CJ FCAPPl M5)*n0.1l 
2Cl tCNUAUI 
C 
C 
C 
C 

C FIND PICOLE OF SYPOCLpCUtR HPNGE 
5VPPPM»P*KSYP 
SYHPl*5YHHP*/2,e0C ♦ C.!UO 
C 

C CCHPUII NCHBER OF HCRUCNIAL LINES |H HUP 
NROSEC>N«PD 
RPCSCC>NRPO 
N«Y«4*NH0SEC/5 

IE INPHCTC .EQ. II PHY>PHOSEC 
NY>?*NRY ♦ I 
NRYl»MLY ♦ I 
R Y I *KR Y 1 

C FINE LEEf AND RIGHT UHII5 OF MAP EY CALLING CIRC 
CALL CIRCINRCSEC) 

C 

DC IC NDMNAP.1,1 

C 

C INItlPLUE LINE VALLES 
C KMAP IS FCR PRINT-CUT PICTURE 
C TRAP IS ECR DICCMEO PEOIC 
CC A K«lil2l 
TPAPUl « RETAim 
A YMAPiRI«ALPHP(23l 
C 

C SPACE CCHN THE PACE 
WRITE 16)161 
It FCRMAI mill 

N5PPCE-I2B - NRY 
OC 17 Klt»l)N5PACE 
17 WRITE I6)IR| ALPHAi2)l 
le FCRPAT I10A,IAU 
C 

C FINE IFE RELATIVE DENSITY, DENSITY, CR DENSITY VARIATION CN EACH 

c hcrucntal line 

DC 2 LP.l,NY 
Ll-AY -LP » I 
FLl«L2 

KPIK»LOWHIILP)ll 

RMAIULOWHIILP,2l 

C 

OC 3 R«RMtN,KMA)l 

C FINE THE K,Y,{ CDCRCIPATES OF EACH POINT ON THE LINE 
FR»R 

liP*IFK - 6l.0C0l*IRLPPER/RADSEC) 

STRECh>5.0DO/A.OOC 
IF (KPHCTO ,E0. II STRECH'I.OCO 
YP«|FL2 - RYII*ISTRfCH*RUPPERI/RADSEC 
R>(CCSOPGI*XP - ICaSXt*SlNDMCI*yp 
V*ISINOPCl**P ♦ ICO!AI*COSOMCI»YP 
2<(SINXII*YP 

C FINE THE SPHERICAL POLAR COORDINATES OF POINT 
R>OSCHTI)(P**2 ♦ YPRA2I 
C IS THE POINT IN THE CCRE7 
IF IB - RCOREI 2S,21,2T 
28 JNUP«22 

C 6 IS USEO FCR CORE 
CU TC 2S 
27 CCNIINUE 

C IS THE POINT BELOW RLCWERT 

IF IR - RLOWERI ICl,lC2,t02 
ICl OIN«C,OCO 
CC TC ICS 


SLI 1P2 
Sll 161 
SLI lAA 
Sll 16S 
SLI UP 
SLI UT 
SLI UP 
SLI U4 

SLI no 

Sll III 
SLI IT2 
SLI IT3 
SLI IIA 
SLI ITS 
SLI IT6 
SLI ITT 
Sll ITS 
SLI IT4 
SLI liO 
SLI m 
SLI 1(2 
SLI 1(1 


$11 ISA 
Sll US 
SLI UP 

SLI UT 
SLI IM 
Sir US 
SLI ISO 
SLI ISl 
SLI IS2 
SLI IS3 
SLI ISA 
SLI US 
SLI U6 


SLI ISB 
SLI ISR 
SLI 2C0 
SLI 201 
SLI 202 
SLI 203 


SLI 20A 
SLI 203 
SLI 206 
SLI 207 
SLI 20a 
SLI 20S 
SLI 210 


SLI ZIt 
SLI 212 
SLI 213 
SLI 21A 
SLI 215 
SLI 216 
SLI 217 


SLI 2U 
SLI 2LS 
SLI 220 
SLI 221 
SLI 222 
SLI 223 
SLI 22A 
SLI 22S 
SLI 226 
SLI 227 
SLI 22B 
SLI 22S 
SLI 230 
SLI 231 
SLI 232 
SLI 231 
SLI 21A 
SLI 215 
SLI 216 
SLI 217 
SLI 21S 
SLI 21S 


» 


] 



52 


r»oor>r»or>rt<*»r>r>r»oonoonnnnonnor»r>or>nr»r»r»oorkooooor»innino rv no o o 


c 

c 


IC2 

CCNI INUI 

SLI 

240 


1I$8YSC*CSCRTI***2 ♦ V**JI 

SU 

241 


IP III 24,25.24 

SLI 

2*2 

2S 

THETAiPICIVE 

SU 

241 


OC 1C 26 

HI 

244 

24 

U*X$CVSC/I 

SU 

2*9 


THEIAyOATaNTU 

SLI 

2*6 


IP Ul 44.26,26 

SU 

247 

44 

THITAtPI A theta 

SLI 

2*1 

16 

CCNTINUE 

SU 

249 


CALL ANClCtKtVtXLPOAl 

SLI 

2S0 


CCHPUTE THE RELATIVE EENSITY* DENSITY, OR OENSITV VARIATICN AT 

POINTSLI 

291 


CALL CCNSTY m*THETA,XLH0A,NDEN»RePOEN,O£NI 

SU 

252 

tC} 

CCNTINUE 

SU 

253 


POT THE RIGHT SYPPOL/HUHEIR IN XPAPIKI/TNAPCRI 

$11 

254 


XS«ICEN/SCALEI a C.TOC 

SLI 

255 


JNUPaSVPBL a rs 

SU 

256 


IP tJNUP >LE. I) }NiP*l 

SU 

251 


IP UNUP .CC. PAXSVHI JNOP>PAXSyP 

SLI 

291 

29 

XP*P|X|»ALPHA|JNUPI 

SLI 

259 


TPAPIXI » DETAIJNLNI 

SU 

260 

1 

CCNT INUE 

SU 

2*1 


c 

C MIhT fh£ tlNE 
HRIlf |«iM) 

14 (SKtUiAU 

WRITE I lO.lOOMTNIRtRl, R«I»I7|) 
100 FCRRRT UZU7I 


SU 2*7 
SU 2*2 
SLI 2*4 
SU 2*S 
SU 26* 


RESET THE LINE 
DC IS K«KH|N,KM*X 
XMPIKI«*IPH*|2)I 
TP6PIKI • BETinil 
S CCN1INUE 
CCNTINUE 


SLI 26T 
SLl 26« 
SU 269 
SU 270 
SLI 271 
SU 272 


CCN1INUE 

RCTLRN 

END 

SUBRCUriNE V6N0ER I PRPiHNAPtWPRPfMXPI 


SLI 27} 
SLI 274 
SLI 275 
VAN 1 


4VANCER* 


VAN 2 


IMS SUMO'JilNE CCNPUei tHL RELATIVE OENSITVi OENSiTY, C« VAN } 

CEnSITY variation on a surface of CHCSEN radius using the van OCR VAN 4 

CRINTEN PROJECTICN. VAN 5 


IT CIVICES THE PAP INTO 

TWO HEHI SPHERES TO PARE THE PAP BIO 

VAN 

6 

ENOUGH 1C OVERLAY CN PALL LCWPAN’S TECTONIC A TIVITY PAP. 

VAN 

7 

NCTATION 



VAN 

e 

RXP • RADIAL OISIANCE OF PAP 

SURFACE IN KILOPETCRS 

VAN 

9 

HPAP • TOTAL HEIGHT Of OVERLAY 

PAP IN CENTIPETERS 

VAN 

10 

WRAP . TOTAL WICTH CF 

OVERLAY 

PAP IN CENTIPETERS 

VAN 

11 

NRAC • RACIUS CF HAP 

IN HCRIICNTAL SPACES ACROSS THE PACE 

VAN 

12 

IPHCTCI 



VAN 

13 

5UBRCUTINES CALLECI 



VAN 

14 

1. CIRC 



VAN 

15 

2. CUBIC 



VAN 

16 

3. CENSTY 



VAN 

17 

4. CUACRC 



VAN 

IB 

nctes 



VAN 

19 

1. NR4C .LE. 102. 



VAN 

20 

2. hPAP .LE. 52 CP. 



VAN 

21 

SAPFLC INPU? DATA ICOlUPN 1 OF 

INPUT STARTS IN CCLLPN 3 HERE) 

VAN 

22 

SUBNPE 


IPORPATT A6I 

VAN 

23 

VANCER 



VAN 

2* 

ICCPENTUI, JAI.141 


IFORPATi l|A6,A2| 

VAN 

25 

Surface pap 



VAN 

26 

NRAD RKP HPAP 

WRAP 

IFORPAH IS.3FI0.5I 

VAN 

27 

102 6371.0 3e.2 

52.0 


VAN 

2« 


IHPLICII REAL*6IA-H)0-n 
INTECERA2 THAPlUt), EEIAIIII 


VAN 29 
VAR JC 


S3 




VAN 

J1 

CI^CKSICN XK*rU2n 


VAN 

12 

CI>«t^SICK CCftMTIMI 


VAN 

}J 

ec>'»ch/eu*/)i»'*f 


VAN 

}A 

CC»'^CK/eUC/f»^l.>'IJ»>ltlVl»5C4U,NPHOrO,M**SyK 


VAN 

J5 



VAN 

It 



VAN 

IT 

CC**»'CN/HKk/liAC*ye,PCCR(*Rt,CWCltfRUPI»Eii 


VAN 

If 

CCRRCMUKI/CCRtM 


VAN 

14 

CrRRCk/5lXJ/kCfA*RffCCN 


VAN 

AC 

cct'»'C^/rur/BllA 


VAN 

At 

ihPir 0*1* 


VAN 

A2 

hRIlf Uf|6) 


VAN 

AJ 

HRtIC >t»]00> 


VAN 

AA 

)t0 FCRktr l///)ICXi'SUfRCU1INE CALUOl VANOER'f///> 


VAN 

AS 

WHITE ItOtll leOHEATUl* Jxl.lAt 


VAN 

At 

)]i rcHHAT iioit»nAe«A2,A//i 


VAN 

AT 

WHITE It. 1011 HHW 


van 

At 

JCI fCHWAI U0*,*HKW ••.lA,rie.3.1H.»HW».10)t»'(HAOlAt BlSIAACe 

OF SURFVAN 

A4 

.ACfl*.//) 


VAN 

SO 

WHITE ltf}02) NHAC 


VAN 

51 

ICJ fCHHAI I10)(.»NHA0 »*.I)(.M.IO)(«MH*OIUS OE WAH |H MOKIUNIAL SHACIVAA 

52 



VAN 

SI 

r<RK|./AAOHVC 


VAN 

5A 

; ElNC HlCCUl Of SVWBCE/CtEM HAWGE 


VAN 

$$ 

SVHf'AHxJ'AASyW 


VAN 

St 

$yHet»SVHWAK/2.C0C « C.lOO 


VAN 

5? 

; CCWHU»e klHHEH CE kCH|<CWIAL UhES IN MAP 


VAN 

St 

NHmCxARAD 


VAN 

59 

KA0$EC>HHAU 


VAN 

to 

NHV»A*NHp5EC/5 


VAN 

61 

IF IMHCrO .EC. II AHyxARDSEC 


VAN 

62 

NyxJAKHV * 1 


VAN 

6) 

NHyi»NRy ♦ I 


VAN 

6A 

RyixARyt 


VAN 

6S 

FIND lEFT ANO RICHT UMIII OF MAP By CAUIMO CIRC 


VAN 

66 

CAEl CIHCINHCSECI 


VAN 

6T 

CC 1C NEMMAP»l.l 


VAN 

6t 

: DC WESTERN HEMISPI'ERC. THIN EASTERN HEMISPMERE 


VAN 

64 

CC 1C NMEMISxl.2 


VAN 

TC 

INITIAUtE LINE VALtEE 


VAN 

Tl 

RMAF IS FCR PRINT-CCT PICTURE 


VAN 

T2 

TRAP IS FCR OICCMED PFOTO 


VAN 

Tl 

DC A K*l.l2t 


VAN 

TA 

tPAPIRI - RETAIIll 


VAN 

T5 

A XMAP(KI»AIPHAI23) 


VAN 

T6 

WRITE It. 161 


VAN 

TT 

It FCRPAT lUUI 


VAN 

Tl 

NV-IRYl - l.CCOl»PM/P/WPAP 


VAN 

T4 

LFMINxNHy - NV 


VAN 

to 

LPMAR-NT - IPPIN * 1 


VAN 

11 

; FINC TEE RELATIVE CfNIITY. DENSITY. OR DENSITY VARIATION CN 

EACH 

VAN 

12 

; HCRIKNTAL LINE 


VAN 

ei 

OC 2 LP«LPMIN,LPMAIf 


VAN 

BA 

IF INFEPIS .EG. 21 CO TC SI 


VAN 

as 

TERpxO.COO 


VAN 

16 

KPINxI 


VAN 

IT 

KMAR>LQWHIILP.2| > <0 


van 

II 

CC TC S2 


VAN 

14 

51 TERPxRACSEC 


VAN 

40 

RPlNxlOWHIILP.il - to ♦ NR05EC 


VAN 

41 

RMAR’NRCSEC 


VAN 

42 

52 CCNTINUE 


VAN 

41 

LI-NY -IP ♦ 1 


VAN 

4A 

FLtxLt 


VAN 

«5 

DC 3 RxRMIN.RMAR 


VAN 

46 

; FIND THE R.V.t CDCRCINATE! OF EACH POINT ON THE LINE 


VAN 

4T 

FRx* 


VAN 

41 

R.FR - l.OOO - TERM 


VAN 

44 

STRECH«5.000/A.ODC 


VAN 

100 

IF INPHCTO .EC. II STHECHxl.DCO 


VAN 

101 

V«(FL2 • RYIIPSTRECF 


VAN 

102 

FINC THE SPHERICAL POLAR COORDINATES OF POINT 


VAN 

10] 

CALL CUeiCIR.Y.RACSIC.TFETAI 


VAN 

lOA 

CALL CUAORCIX.Y.RAOSEC.ALPOAI 


VAN 

lOS 

CCMPUTE THE RELATIVE CENSITV. OENSITY. OR CENSITY VARlATICN 

AT POINTVAN 

106 

CALL CENSTV IR. THETA. aLPOA.NOEN.REFOEN.OENl 


VAN 

lOT 

PUT THE RIGHT SVPBOL/FUPBER IN RMAPIR l/TMAPI K 1 


VAN 

101 

liSxICENASCALEI • C.'OC 


VAN 

104 

JNUPxSVPRL « AS 


VAN 

no 


54 



IF tJNUR .LE. 1) JNIH>1 



VAN 

111 


ir UftUR .CE. HAX5VM JAUH«HA»SyR 



VAN 

112 

29 

RRAP4K|i>ALPHA<JNUR> 



VAN 

113 


TRAPltU « BETAIJNIRI 



VAN 

UA 

) 

CCNTIRUE 



VAN 

115 

C 

PRIRT THE LINE 



VAN 

IIA 


HRITE giifUl (XPAPIPI. P»t.l2l) 



VAN 

117 

M 

PCRPAT ISX.I2IAI) 



VAN 

HE 

C 

MRITE 1 IDtlCCUTRAPlRl, K«1.I21) 



VAN 

IIV 

C 100 FcnfAT I121I2I 



VAN 

120 

w 

c 

reset the lire 



VAN 

121 


DC 15 K»tiHIN,KPAK 



VAN 

122 


XPAP|R)>ALPHAI23I 



VAN 

123 


THAPIKI • RETAItn 



VAN 

12A 

15 

CCNTIRUE 



VAN 

125 

1# 

2 

#► 

CCNTlRUe 



VAN 

126 

I# 

1C 

CCNTIRO^ 



VAN 

127 


RETIRR 



VAN 

12B 


ERD 



VAN 

129 

C 

SUeRCUiMRC OERSTY 1 F| THETA, XLHOA,NDEN|REFDEN.OENI 


DEN 

1 

c 

c 

THIS 5UER0UTIRE CCRfUIES THE RELATIVE OERSITY, 

DENSITY, CR 

DEN 

2 

C CENSnV VtRUTlON AT * fOINT INSICE THE EARTH. 
C 



DEN 

3 

C 

c 

r 

RCIATIOR 



DEN 

A 

V 

c 

R • RACIAL 0IS1ARCE OF POINT FROH THE CENTER 

CF THE EARTH, 

CEN 

5 

c 

N0RPALI2ED TC THE EARTH'S RADIUS 



DEN 

6 

c 

THETA « COLATITUCE OF FOINT IN RAOIANS 



DEN 

T 

c 

XLPCA • ANCLE CF PCIRT WITH THE X-AXIS, LYING 

IN 

THE X-Y PLANE, 

DEN 

B 

c 

HEASURED IR RAEIANS 



DEN 

9 

c 

DIN • relative density, DENSITY, OR DENSITY 

VARIATION AT THE 

DEN 

10 

c 

POINT 



DEN 

11 

c 

NCER « 0 FOR RELATIVE DENSITY, >1 FOR ACTUAL 

DENSITY, 

DEN 

12 

c 

c 

-2 FCR DENSITY VARIATION 



CEN 

13 

c 

c 

e 

SUBRCUTINES CALLED! 



DEN 

lA 

c 

1. LCERDR 



DEN 

15 

c 

c 

2. C2IEUN IIF RDEN'l OR N0EN>2I 



DEN 

16 

c 

IPPLICIT flEAL«6(A-H,0-n 



DEN 

17 


CIPERSICN CLHI36,3T,2l,COE|36) 



DEN 

IB 


OIPERSICN PBARIT03) 



DEN 

19 


CCPHCR/CLKH/CLP.CCE 



DEN 

20 


CCFFCR/ELKC/F,PI,Pi:,PltlV2,SCALE,NPHOTO,MAXSYH 


CEN 

21 


CCHHCN/eLKCALFIR,LNAX,LlOH,LUP 



CEN 

22 


CCHHCN/RLRH/RAOAVG,PCCRE,RLCMER RUPPER 



DEN 

23 


CCNPCR/eLKF/PaAR 



DEN 

2A 

c 

CHANCE VARIABLES FRCN THETA TC PHI 



DEN 

25 


PHI>PI0IV2 - THETA 



DEN 

26 

c 

CCMPUTE LEGENDRE POLYROPIALS 



CEN 

27 


CALL LCENDR(PHI,LPAX) 



DEN 

28 

c 

CCNPUTE RELATIVE DENSITY 



DEN 

29 


SUK>C.OCO 



CEN 

3C 


DC 5 L»LPIN,LPAX 



DEN 

31 


RL-R»»L 



DEN 

32 


MAX.L ♦ 1 



DEN 

33 


DC 5 PI«1,NAX 



DEN 

3A 


XH*P1 - 1 



DEN 

35 


arg<xp*xlmda 



DEN 

36 


SS'CSINIARGI 



DEN 

37 


CCaCCCSIARC) 



DEN 

38 


IND<P1 * L*(L*ll/2 



DEN 

39 


SUH»SUH F ICLRIL,P1,U*CC ' CLHIL,NI,2)*SSI*IPBARI INO) 1 *RL 

DEN 

AO 

5 

CCNTINUE 



DEN 

Al 


OEN«SUN 



DEN 

A2 

c 

SEE IF ACTUAL CENSITY IS HANTEO 



DEN 

A3 


IF INCEN .EC. 1) CO TC T 



DEN 

AA 

c 

SEE IF DENSITY VARIATION IS WANTED 



DEN 

AS 


IF INCEN .Efl. 21 CO TC S 



DEN 

A6 


CC TC B 



DEN 

A7 

f 

CCNTINUE 



DEN 

AB 


CALL CIIEHN IR,OENOIII 



DEN 

A9 


0EN>CENC2I*I l.OCO * DEN 1 - REFDEN 



DEN 

SO 


CC TC S 



CEN 

SI 

9 

CCNTINUE 



DEN 

52 


CALL C2IEHN IR,DENOm 



DEN 

S3 


DEN>CENC2I*0EN 



DEN 

SA 

6 

CCNTINUE 



DEN 

55 


RETURN 



DEN 

56 


55 


r.%u 

SUBKCiiUMB 


nih s? 
m I 


C TMi SUBAtiunKC CCNtUtI{ tHB ANCLE WIIM tllC X'AfSIS EMON DIE lit«VI 
C CCCWCINAUS OE * ECINJ. 

c 

c 

C EC1AIIDN 

C 

e * • X-CCOACIAAIC Of EOlNt 

c V • y-cccM5iK*ie of coini 

C XLNCA • ANCLE NtfH A>AA|$ 


SUNNEOflNES CALLECi 


t. ALNDA IS NEASUNEL IN AADIANS. 
2. tr X<Q AND V«Qf DIEN KLHOA«Q. 


lAELICM NEAL*«|A-ll,0-EI 

H»E.UI»‘f!*CS5EyfStO 

PIOIV3«EI/E.OCO 

If IK) >1.33.21 

U«¥3* 

XLACA<t:ArANIUI 
If l¥) 3A«3St3A 
ALHCAaXLNOA « IM 
CC 1C 10 
XLNCA»PI 
GC 1C 30 
CCNMNllE 
If IV) 27.3H.3<I 
NLNCA— flUIVl 
OC 1C 10 
XLfCA*O.OnO 
ee TC 30 
XINCA*I'IDIV3 
CC K 10 
CCNIINUE 
If lYI 11.13.11 

u*v/* 

XLNCA'CATANIUI 

CC 1C 30 

XLNCA'O.OOO 

CCNIINUE 

HETUAN 

END 

SUHACUIING LCENCHtPM.NfAX) 


C THIS SUHHOUIINE CCNI'UTEJ THE NDKfAUlEO ASSOCIATED LECENDHE 
C POUNCHIALS. 

C 

c 

C NCTATION 

t 

C PHI • CCOCENTHIC LAI |TUOE IN AAOIANS 

C NPAX • NAXIPUP CCCHEC AND OHOEH Of THE MOLTNOHIALS 

C PPAP • ANAAV IN NfICH DIE POLTNOHIALS AHC SICKED 


SUAACUTINES CALLECI 


1. CCHAUTES 1‘CLVNOflALl USING PAULA'S A»PI NQAPALI2ATICN. 

2. NPAX IS LESS THAN OP EClAL TC IAO. 

S. CIPINSION Cf PEAR IS I ♦ NPA*»|NPA«»ll/3. 

A, PCLYNOnlAL CP UICRCE N AND ORDER P IS STORED IN PnARlINOeiO. 

NHBR6 INCeX . I ♦ P ♦ N*INMI/2. 

S. SUSROUTINE HAS LNCERPLONS PHEN PHI IS NEAR IHE POLES. 


IPPLICIT REAL** lA-f.C-Jl 
LCCtCAL NOT 1ST 

DIPENSICN PR((Rnoil,SCRCOm6U 

CCPPCN/ELRf/f'RAR 

DATA NOnSrA.fALSI.A 


ANQ 

2 

ANC 

3 

ANC 

A 

ANC 

S 

ANC 

A 

ANC 

T 

ANC 

E 

ANC 

S 

ANC 

10 

ANC 

\\ 

ANC 

12 

ANC 

11 

ANC 

lA 

ANC 

15 

ANC 

t« 

ANC 

U 

ANC 

IP 

ANC 

IN 

ANC 

20 

ANC 

21 

ANC 

22 

ANC 

21 

ANC 

3A 

ANC 

25 

ANC 

36 

ANC 

2T 

ANC 

26 

ANC 

2N 

ANC 

30 

ANG 

11 

ANO 

12 

ANC 

11 

ANC 

lA 

ANC 

15 

ANG 

36 

ANC 

IT 

ANG 

IB 

ANC 

IN 

ANC 

AO 

LCE 

1 

LOE 

2 

LCE 

1 

LCE 

A 

LCE 

5 

LOU 

6 

LCE 

1 

LCE 

R 

LCE 

N 

LCE 

to 

LCE 

11 

LCE 

12 

LCE 

11 

LCE 

lA 

LCE 

IS 

LCE 

16 

LCE 

17 

LCE 

16 

LCE 

IN 

LCE 

20 

LOG 

21 


8 « 


U' <Hrmi CO 1C I 
Nnm«ii»vc* 

KSUI'*J»N«*K*l 
n (KStCO.Q1.lfcl) St(?P f 

cc * i*uNS(r)> 

fc sst»teiu».oscHi(OH.tfci) in 

fc SIHOM-CSINU'HU 

ccsPM»ccusm'ii 

OOlMINCKOlQiOiniOC 

oB*MiNuiu'n.Qu*SQPeci m»siM'ni 
OBAOi INCXOU. nitSOPQCT (ii*co$oin 
CC ?c 

PRAOnNCXlMN.CiniQPCCI l!*Nn)/C»lO*T(Nl« 

. (SlK(On::*N«U*SINOHt*OOfcMINO)(0(N«l.Dll 

, -lUUfcltN-VlfcORlMINOKlMN-JltOn/SCOCOnJloX-in 

u«t5coi;i;tii*K»n*com I 
u«stRccm)«ii*i*r»M iNcxotN-nen 

If ll-NI 

B tl»U*SQHOOnK*n*SUt01tl»-?l»f'H»tmNa)(IMN-?»Ul«CBOOI (J»N’1I 
>l OBfcMINCXNN.UnsePCCI (?»N*U»U/1S«BCO(INI»seBCOIIN*n) 
nt* ?c o.j.N 

ij»u«oB*BnMiKPiK-i.o-in 
If (f-N*3!l 11(11.17 

n U-l7*seBClnifc-l*l*SCRCOUN-H-n»OHfcM INt))l0lN»7,Hn/SCBCCtl7*fcl“ll 
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